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SUMMARY 


The  use  of  ceramics  to  construct  underwater  pres¬ 
sure  resistant  housings  requires  technique 
can  be  used  to  adjoin  ceramic  hull  compc 
together.  The  Naval  Command,  Control  a 
Surveillance  Center  (NCCOSC)  RDT&E  Div.__,. 
(NRaD)  has  pioneered  the  use  of  epoxy-bonded 
metallic  joint  rings  for  the  purpose  of  assembling 
adjacent  ceramic  housing  sections.  These  joint 
rings  act  to  transfer  load  through  joint  interfaces, 
while  also  providing  a  means  of  sealing  and  main¬ 
taining  closure.  Additionally,  joint  rings  can  be  de¬ 
signed  to  provide  additional  stiffness  to  the  housing 


assembly  to  increase  buckling  resistance. 

This  report  summarizes  issues  that  must  be  con¬ 
sidered  when  designing  metallic  joint  rings  for 
ceramic  underwater  housings.  Selection  of  the  joint 
ring  material  and  joint  ring  bonding  adhesives  are 
addressed.  The  effect  of  various  joint  ring  design 
parameters  on  the  structural  performance  of  the 
ceramic  housing  assembly  are  also  discussed.  The 
report  concludes  with  a  description  of  four  different 
alumina-ceramic  pressure  housings  that  were 
designed  and  pressure  tested  to  destruction.  The 
purpose  of  these  tests  was  to  validate  the  structur¬ 
al  performance  of  various  joint  ring  designs  for 
consideration  in  future  ceramic  underwater  pres¬ 
sure  resistant  housings. 
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INTRODUCTION 


Figure  1 1  shows  three  cross  sections  of  designs 
that  could  be  used  to  assemble  a  ceramic  cylindri¬ 
cal  hull  for  an  underwater  pressure  housing.  These 
joining  techniques  differ  in  the  configuration  of 
metallic  joint  rings  that  are  used  at  the  interface  of 
adjacent  ceramic  components.  Joint  rings  may 
serve  a  variety  of  different  functions  to  the  hous¬ 
ing,  but  all  share  certain  primary  functions  essen¬ 
tial  to  their  successful  application  in  ceramic 
underwater  pressure  housings. 

Joint  rings  act  to  transfer  load  between  adjacent 
ceramic  sections  without  having  the  bearing  sur¬ 
faces  at  the  end  of  one  ceramic  part  bear  directly 
against  the  bearing  surface  of  an  adjacent  ceramic 
part.  Early  attempts  to  join  ceramic  components 
without  the  protection  of  joint  rings  (reference  11) 
resulted  in  short-term  cyclic  fatigue  failure  initiated 
by  cracks  that  originated  from  the  region  of  direct 
ceramic-to-ceramic  contact.  Protecting  the  ends  of 
ceramic  components  with  metallic  joint  rings  allows 
transfer  of  pressure-induced  interface  loads 
between  hull  sections  without  causing  fretting  fail¬ 
ure  of  ceramic  ends. 

Joint  rings  also  provide  local  attachment  points  for 
mounting  internal  or  external  hardware  such  as  tie 
rods,  payload  rails,  and  electrical  cable  raceways. 
The  means  for  achieving  joint  seals  and  joint  clo¬ 
sure  can  also  be  designed  into  joint  rings.  By  inte¬ 
grating  these  features  into  the  joint  ring,  the 
localized  loads  associated  with  handling,  assem¬ 
bling,  and  sealing  the  pressure  housing  can  be 
directed  onto  the  metallic  end  rings  and  not  on  the 
more  sensitive  ceramic  hull  components. 

Additionally,  it  may  be  advantageous  to  design  the 
joint  rings  to  provide  additional  stiffness  to  increase 
the  buckling  resistance  of  the  pressure-housing 
assembly.  Manufacturing  constraints  may  limit 
techniques  that  the  designer  can  use  to  stiffen  indi¬ 
vidual  ceramic  hull  components.  In  these  cases, 
metallic  joint  rings  can  be  configured  to  provide  the 


1 .  Figures  and  tables  are  placed  at  the  end  of  the 
text. 


required  ceramic  component  end  support  needed 
to  achieve  the  operating  depths  for  which  the 
housing  is  intended. 

JOINT  RING  FUNCTION  _ 


The  first  housing  configuration  shown  in  figure  1 
would  appear  to  be  the  simplest  way  to  construct  a 
ceramic  cylindrical  hull  section  for  an  underwater 
pressure  housing.  This  assembly  consists  of  a 
single  monocoque  ceramic  cylinder  with  appropri¬ 
ate  metallic  end  caps  for  mating  with  the  end  clo¬ 
sures.  But  for  many  housing  configurations,  a 
single  monocoque  cylinder  hull  may  not  be  a  viable 
design  approach.  Large  housings  (reference  6), 
with  outer  diameters  of  25  inches  or  larger  and  hull 
lengths  corresponding  to  L/D  ratios  greater  than 
two,  may  simply  not  be  manufacturable.  Fabrica¬ 
tion  of  these  large  parts  may  be  limited  by  the  size 
of  available  isopress  equipment  kilns,  and  grinding 
machines.  Additionally,  cylinder  length  may  be 
limited  by  other  constraints  such  as  slumping  of 
the  green  body  during  firing.  Successful  construc¬ 
tion  of  such  large  cylindrical  housings  requires 
assembling  a  number  of  shorter  cylinders  together 
using  joint  rings  or  other  bonding  techniques. 

Brazing  a  number  of  ceramic  rings  together  has 
been  demonstrated  to  be  one  potential  approach 
to  building  up  longer  cylindrical  hull  sections  that 
may  not  be  fabricated  as  one  piece  (reference  11). 
Six-inch  outer-diameter  (OD)  cylinders  (L/D  =  1 .5) 
have  been  constructed  by  metallizing  the  ends  of 
ceramic  rings  and  then  brazing  the  rings  together. 
The  shortcoming  of  this  and  similar  techniques 
such  as  diffusion  bonding  is  that  the  size  of  cylin¬ 
der  that  can  be  constructed  this  way  is  limited  by 
the  size  of  the  fabrication  equipment  that  is  avail¬ 
able.  Brazing  furnaces  that  are  big  enough  to  han¬ 
dle  larger  ceramic  components  may  not  exist  or  be 
of  limited  access. 

Another  technique  for  joining  ceramic  cylinders 
involves  using  a  metallic  coupling  ring  like  that 
shown  in  the  second  housing  configuration  in  fig¬ 
ure  1 .  This  technique  involves  epoxy  bonding  an 
H-shaped  ring  to  the  ends  of  adjacent  ceramic  cyl¬ 
inders.  Assembling  cylinders  in  this  way  can  be 
performed  using  relatively  simple  handling  fixtures 
and  is  not  size  limited.  Pressure  testing  of 
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cylindrical  hulls  constructed  using  coupling  rings 
has  demonstrated  the  viability  of  this  approach. 
Analysis  and  testing  has  shown  that  the  structural 
performance  of  the  coupling  ring  configuration 
matches  or  exceeds  the  performance  of  the  single 
monocoque  cylinder  it  replaces.  As  mentioned  pre¬ 
viously,  the  presence  of  coupling  rings  also  pro¬ 
vides  a  means  of  accommodating  local  internal 
and  external  attachments  to  the  pressure  housing 
walls.  The  general  structural  behavior  of  ceramics 
is  such  that  they  are  well  suited  to  bearing  the  pri¬ 
marily  compressive  membrane  stresses  that  exist 
in  the  walls  of  underwater  pressure  housings  under 
depth  loading.  On  the  other  hand,  ceramics  will  not 
perform  as  well  in  areas  where  localized  stresses 
exist  due  to  local  attachments  and/or  joint  inter¬ 
faces.  For  this  reason,  metallic  coupling  rings  are 
advantageous  in  helping  to  buffer  the  ceramic  hull 
against  these  extraneous  stresses. 

In  addition  to  manufacturing  considerations  that 
may  limit  their  use,  there  are  a  number  of  housing 
configurations  where  a  monocoque  cylinder  design 
is  undesirable  for  structural  reasons.  The  structural 
integrity  of  cylinders  with  large  values  of  L/D  sub¬ 
jected  to  high  hydrostatic  pressure  may  depend  on 
their  ability  to  resist  buckling.  For  this  type  of 
monocoque  hull,  the  wall  thickness  that  is  required 
to  resist  buckling  may  have  to  be  substantial. 

While  the  thick  hull  may  be  capable  of  avoiding 
collapse  by  buckling,  it  could  be  understressed  to 
the  point  that  the  high  compressive  strength  of  the 
ceramic  material  is  not  being  utilized.  This 
approach  would  result  in  a  housing  with  a  higher 
weight-to-displacement  (W/D)  ratio  and  would  give 
lower  performance  than  could  be  achieved  with 
other  housing  designs. 

Designing  a  more  efficient  housing  by  using  inte¬ 
gral  stiffening  ceramic  ribs  would  be  a  more  effi¬ 
cient  approach  from  a  weight-savings  point  of  view, 
but  unfortunately  is  restricted  because  fabrication 
is  not  economical  because  of  manufacturing 
constraints.  To  have  stiffening  ribs  in  large  cylin¬ 
ders  made  from  ceramic  materials  such  as  alu¬ 
mina,  one  would  have  to  isopress  a  very  thick 
walled  cylinder,  and  then  create  ribs  by  grinding 
out  the  excess  material  between  the  ribs  once  the 
ceramic  body  has  been  fired.  Th,s  would  be  an 


exorbitantly  expensive  and  risky  operation.  Green 
machining  this  excess  material  prior  to  firing  is  not 
an  option  because  of  thermal-expansion  mismatch 
problems  and  the  potential  for  high  residual 
stresses  that  would  occur  at  the  rib/hull  interface 
during  firing.  A  typical  alumina-ceramic  green  body 
that  is  fabricated  by  isopressing  can  shrink  as 
much  as  20%  during  the  firing  process. 

The  prospect  for  integral  ribs  is  much  better  for 
cermet  hull  components  made  by  alternative  fab¬ 
rication  techniques  (reference  12).  An  example  of 
this  would  be  silicon  carbide-reinforced  alumina- 
ceramic  matrix  composites  made  by  directed  metal 
oxidation  process.  This  material  process  results  in 
less  than  a  one-percent  dimensional  change  during 
fabrication  which,  in  theory,  should  allow  for  near 
net-shape  fabrication  of  housing  components  with 
features  like  integral  ribs.  An  additional  benefit  that 
cermet  materials  could  have  for  underwater  pres¬ 
sure  housings  is  that  the  shell  wall  could  transition 
to  a  metal-rich  composition  toward  the  ends  of  the 
hull  section.  This  would  allow  the  cermet  cylinder 
ends  to  be  designed  with  features  such  as  O-ring 
glands,  attachment  points,  and  flanges  for  closure. 
The  size  and  shape  of  housings  made  to  date  by 
this  process  are  limited  to  12-inch  OD  monocoque 
cylinders  with  L/D  =  1 .5. 

Given  these  current  limitations  in  fabricating  stiff¬ 
ened  ceramic  hull  components,  alternative 
approaches  to  increasing  the  buckling  resistance 
of  ceramic  housings  are  required.  One  such  tech¬ 
nique  is  to  integrate  stiff  metallic  joint  rings  into  the 
housing  assembly  design.  These  rings  can  be 
designed  to  support  the  cylinder  ends  in  such  a 
way  as  to  raise  the  collapse  pressure  of  the  hous¬ 
ing  by  increasing  the  external  pressure  required  to 
bend  the  ceramic  hull  into  a  buckled  configuration. 
The  last  of  the  three  housing  configurations  shown 
in  figure  1  uses  a  central  metallic  joint  ring  assem¬ 
bly  designed  to  stiffen  the  ceramic  hulls.  The  use 
of  stiff  joint  rings  to  help  bear  the  external  hydro¬ 
static  load  allows  the  designer  to  build  up  cylindri¬ 
cal  hulls  by  using  a  number  of  cylinders  which 
have  a  reduced  wall  thickness  compared  to  the 
single  monocoque  cylinder  they  replace.  This  com¬ 
posite  structure  of  stiff  metallic  joint  rings  and 
ceramic  hulls  allows  the  designer  to  achieve  the 
required  buckling  resistance  of  the  assembly  while 
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tailoring  down  the  ceramic  wall  thickness  to  utilize 
ceramic’s  excellent  compressive  strength.  This 
allows  structures  to  be  optimized  to  have  a  mini¬ 
mum  W/D  ratio  by  taking  advantage  of  the  material 
properties  of  each  of  the  housing’s  components. 

JOINT  RING  MATERIAL 


In  the  design  of  joint  rings  for  underwater  ceramic 
pressure  housings,  material  selection  is  driven  by 
the  specific  criteria  that  govern  the  individual  hous¬ 
ing  design.  The  use  of  ceramics  is  attractive  for 
underwater  housing  designs  where  low  W/D  ratios 
are  needed  to  provide  maximum  buoyancy. 
Ceramics  are  an  excellent  candidate  for  this  pur¬ 
pose  because  of  their  high  specific  compressive 
strength  and  specific  modulus.  Consequently,  the 
selection  of  a  joint  ring  material  for  these  types  of 
applications  would  require  metals  that  also  have 
high  specific  strengths  to  help  keep  the  structural 
weight  of  the  housing  assembly  to  a  minimum.  In 
cases  where  the  metallic  joint  rings  are  designed 
to  provide  additional  stiffness  to  increase  the  buck¬ 
ling  resistance  of  the  housing  assembly,  the  spe¬ 
cific  modulus  of  the  joint  ring  material  becomes  an 
important  concern. 

An  additional  benefit  of  using  ceramic  for  major 
underwater  hull  components  is  its  excellent  resis¬ 
tance  to  corrosion  in  seawater.  Likewise,  materials 
used  for  joint  rings  should  be  capable  of  withstand¬ 
ing  the  marine  environment  without  degradation 
due  to  any  type  of  corrosive  failure. 

A  number  of  ceramic  materials  like  beryllia  exhibit 
outstanding  heat  conductivity  which  make  them 
attractive  for  applications  in  housings  that  require 
heat  dissipation  to  keep  internally  packaged  com¬ 
ponents  such  as  electronics  at  low  ambient  tem¬ 
peratures.  For  these  types  of  applications,  it  may 
also  be  desirable  to  select  materials  with  high  ther¬ 
mal  conductivities  to  fabricate  the  joint  rings  in 
order  to  offer  an  additional  path  through  which 
internally  generated  heat  can  be  dissipated. 

In  situations  where  underwater  housings  require 
the  capability  to  operate  over  a  wide  range  of  tem¬ 
peratures,  matching  the  thermal  expansion  coeffi¬ 
cients  between  the  ceramic  used  in  the  hull  and 
the  material  used  in  the  joint  ring  becomes  impor¬ 


tant.  The  closer  the  coefficents  are  to  each  other, 
the  less  thermally  induced  stress  will  occur  at  the 
interface  between  the  two  materials  under  thermal 
cycling.  Stresses  induced  by  thermal  loading  may 
be  of  concern  specifically  in  cases  where  the  bear¬ 
ing  surfaces  of  the  ceramic  ends  are  placed  in 
localized  tension.  In  such  cases,  the  potential  for 
cracks  to  propagate  from  any  local  tensile  regions 
would  exist  if  the  housing  was  subjected  to  high 
numbers  of  thermal  cycles.  If  high  temperatures 
are  part  of  the  operating  environment  for  the 
ceramic  housing,  special  care  must  be  given  to  the 
type  of  materials  used  for  the  joint  rings  as  well  as 
the  type  of  bonding  adhesives  and  seals  that  are 
selected. 

Another  localized  stress  that  may  appear  at  the 
interface  between  the  ceramic  hull  and  metal  joint 
ring  is  caused  by  a  Poisson’s-type  effect  that 
occurs  under  hydrostatic  pressure  loading  of  the 
housing  assembly  (reference  7).  External  pressure 
on  underwater  housings  results  in  substantial  com¬ 
pressive  membrane  strains  in  the  housing  compo¬ 
nents.  This  meridional  and  circumferential 
compression  leads  to  radial  expansion  of  the  hous¬ 
ing  wall  under  depth  load.  The  amount  of  radial 
expansion  depends  on  the  elastic  moduli  and  Pois- 
sion’s  ratio  of  the  housing  material  and  the  magni¬ 
tude  of  the  membrane  stresses.  At  discontinuities 
in  the  housing-wall  material  and  geometry  such  as 
at  the  interface  between  a  ceramic  hull  and  a 
metallic  end  cap,  a  mismatch  in  radial  expansion 
will  occur.  Typically,  the  metallic  end  cap  will 
undergo  greater  radial  expansion  than  the  rela¬ 
tively  stiff  adjacent  ceramic  shell  and  thereby  place 
the  ceramic  bearing  surface  region  into  localized 
tension.  This  effect  can  be  of  concern  when  the 
underwater  housing  assembly  is  required  to  com¬ 
plete  a  high  number  of  dive  cycles.  The  tensile 
stresses  present  at  the  ceramic  bearing  surface 
may  be  large  enough  to  initiate  cracks  that  propa¬ 
gate  into  the  ceramic  wall  from  bearing  surface 
flaws  under  repeated  loading. 

In  addition  to  material  properties,  consideration  of 
relative  material  and  manufacturing  costs  for  joint 
rings  for  ceramic  underwater  pressure  housings 
are  obviously  important.  The  machinability  rating 
of  the  material  becomes  an  issue  especially  when 
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fabricating  joint  ringc  ,^r  large  housings  (refer¬ 
ence  6).  Ceramic  housing  joint  rings  used  at  the 
Naval  Command,  Control  and  Ocean  Surveillance 
Center  (NCCOSC)  RDT&E  Division  (NRaD)  are 
machined  from  solid  rolled-ring  forgings  which 
.  *=>  .juire  substantial  material  removal  to  create  the 
U-shaped  joint  ring  used  to  encapsulate  the 
ceramic  ends.  In  *oint  ring  designs  requiring  exter¬ 
nal  or  internal  lugs  or  ears  for  attachment  points, 
the  weldability  of  the  material  may  become  an 
important  criteria.  Since  the  structural  performance 
of  pressure  housings  can  be  sensitive  to  geometric 
imperfections,  welding  techniques  that  minimize 
distortion  in  the  finished  parts  should  be  pursued. 

The  choice  of  material  for  use  in  a  joint  ring 
depends  on  the  performance  criteria  selected  for 
each  specific  application.  Both  titanium  alloy 
Ti-6AI-4V  and  high-strength  7000-series  aluminum 
alloys  have  been  successfully  employed  at  NRaD 
for  use  in  joint  rings  for  ceramic  underwater  pres¬ 
sure  housing  assemblies.  Titanium  alloys  are 
especially  attractive  because  of  their  high  specific 
strength,  high  specific  modulus,  corrosion  resis¬ 
tance,  thermal  expansion  coefficient  match  with 
alumina  ceramic,  and  weldability. 

JOINT  RING  BONDING  MATERIALS 

As  mentioned  previously,  there  are  a  number  of 
techniques  that  can  be  used  to  join  adjacent 
ceramic  housing  components.  When  metallic  joint 
rings  are  selected  as  **e  means  of  attachment, 
procedures  for  bonding  the  joint  ring  to  the  ends  of 
each  ceramic  hull  component  must  be  developed. 
Two-part  epoxy  adhesives  had  been  selected  by 
NRaD  for  this  purpose,  and  a  description  of  the 
bonding  procedure  currently  used  by  NRaD  is 
given  in  appendix  A.  The  use  of  two-part  epoxies 
for  joint  ring  assembly  may  not  be  an  acceptable 
choice  in  cases  where  there  are  high  operating 
temperatures.  The  top  assembly  shown  in  figure  2 
shows  a  cross-section  view  of  a  ceramic  cylinder, 
metallic  joint  ring  for  mating  and  sealing  with  an 
end  closure,  and  the  epoxy  bond  filling  the  annular 
and  axial  spaces  between  these  two  parts.  A  gen¬ 
eral  discussion  on  important  considerations  to  be 
used  for  selection  of  a  suitable  bonding  adhesive  is 
as  follows. 


The  bonding  material  should  offer  compliance  to 
eliminate  point  loadings  on  the  bearing  surface  of 
the  ceramic  ends  caused  by  imperfections  in  the 
ceramic  surface.  The  bonding  material  must  be 
capable  of  transferring  the  high  compressive  loads 
between  the  ceramic  hull  and  its  metallic  joint  ring 
without  degradation.  The  joint  ring  design  and  the 
adhesive  selected  must  perform  this  load  transfer 
without  the  adhesive  extruding,  fracturing,  delami¬ 
nating,  or  exhibiting  any  other  type  of  permanent 
failure.  The  use  of  small  chamfers  or  radii  at  the 
edges  of  the  ceramic  hull  component  ends  aids  in 
reducing  the  stresses  on  the  bonding  material  as 
well  as  reducing  the  chance  for  chipping  the  edge 
of  the  ceramic  part  during  assembly. 

The  bond  between  the  joint  ring  and  hull  compo¬ 
nent  must  also  withstand  external  pressure  and  the 
marine  environment  without  any  leakage.  As  dis¬ 
cussed  in  appendix  A,  one  way  of  ensuring  that  no 
water  penetrates  the  epoxy  bond  in  the  joint  is  by 
applying  an  R7V  sealant  over  the  epoxy  where  it  is 
exposed  to  water  as  shown  in  the  top  assembly  of 
figure  2.  This  design  may  be  of  use  when  the  pres¬ 
sure  housing  is  intended  to  undergo  a  large  num¬ 
ber  of  dive  cycles  or  extended  submersion  and  the 
potential  for  water  intrusion  becomes  more  likely. 
Other  techniques  for  sealing  the  epoxy  include 
replacing  a  bead  of  RTV  with  a  urethane  coating  or 
mechanical  seal  like  the  one  shown  in  the  bottom 
assembly  of  figure  2.  This  seal  is  an  elastomeric 
boot  that  is  “rubber  banded”  into  place  and  offers 
the  benefit  of  being  removed  for  easier  inspection 
of  the  underlying  ceramic  hull. 

Another  advantage  already  noted  is  that  joining 
ceramic  components  with  epoxy  bonded  metallic 
joint  rings  is  an  economically  feasible  approach 
that  can  be  performed  using  simple  assembly  fix¬ 
tures  and  is  not  constrained  by  the  size  of  the 
housing  components.  The  existence  of  these  sim¬ 
ple  joining  techniques  that  perform  reliably  during 
service  are  essential  to  the  continued  success  of 
ceramic  underwater  pressure  housings. 

Another  potential  bonding  technique  involves  plac¬ 
ing  a  brazing  alloy  between  the  ceramic  hull  and 
metallic  joint  ring  and  heating  this  assembly  in  a 
furnace  to  obtain  a  brazed  joint.  A  number  of 
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techniques  have  been  developed  for  ceramic-to- 
metai  brazed  joints  including  sintered  metal  pow¬ 
der  process,  active  filler  metal  process,  and  vapor 
coating  process  (reference  9).  A  typical  active  filler 
metal  process  brazing  alloy  requires  a  minimum 
furnace  temperature  of  600  degrees  Celsius  to 
form,  which  requires  a  good  match  in  thermal 
expansion  coefficients  between  the  ceramic  and 
metal  joint  ring  to  avoid  residual  stress  problems 
once  the  brazed  joint  has  been  cooled. 

Residual  stresses  that  result  from  brazing  can 
have  catastrophic  effects  if  their  result  is  to  place 
the  ceramic  hull  ends  in  a  state  of  tensile  stress. 
Failure  by  static  fatigue  due  to  residual  tensile 
stresses  in  the  ceramic  can  occur.  The  larger  the 
diameter  of  the  ceramic  hull  component,  the  larger 
the  potential  for  thermal  deflection  differences  dur¬ 
ing  brazing  and,  consequently,  the  greater  the 
chance  for  residual  stress  problems.  One 
approach  to  dealing  with  thermal  expansion  mis¬ 
match  between  the  ceramic  hull  and  metallic  joint 
ring  is  to  add  an  interlayer  of  a  metal  that  has  an 
intermediate  thermal  expansion  coefficient,  but  this 
necessitates  an  additional  joint  as  well  as  addi¬ 
tional  risk  and  cost.  Other  ceramic-to-metal  bond¬ 
ing  techniques  such  as  fusion  welding,  diffusion 
bonding,  and  glass  sealing  exist,  but  are,  as  yet, 
untested  for  use  in  ceramic  underwater  pressure 
housings. 

JOINT  RING  DESIGN _ 


The  detailed  design  of  metallic  joint  rings  for 
ceramic  underwater  pressure  housings  is  obviously 
dependent  on  the  specific  application  for  which  the 
housing  is  intended.  Nonetheless,  general  recom¬ 
mendations  can  still  be  made  about  certain  details 
the  engineer  will  have  to  consider  when  designing 
metallic  joint  rings. 

Figure  3  shows  three  options  for  joint  rings  based 
on  a  Naval  Ocean  Systems  Center  (NOSC)2  type 


2.  NOSC  is  now  Naval  Command,  Control  and  Ocean 
Surveillance  Center  (NCCOSC)  RDT&E  Division 
(NRaD). 


Mod  1  (hereafter  called  Mod  1)  design  for  a  service 
joint  between  a  ceramic  cylindrical  hull  and  a 
ceramic  hemispherical  end  closure.  The  primary 
functions  of  the  metallic  joint  rings  for  this  case  are 
to  protect  the  bearing  surfaces  of  the  ceramic  ends 
and  to  provide  a  means  of  assembly  and  sealing. 
Protection  of  ceramic  ends  is  especially  important 
when  there  are  substantial  differences  in  radial 
deflection  between  adjacent  hull  components 
under  external  hydrostatic  pressure.  Fretting 
between  unprotected  ceramic  ends  would  be  more 
severe,  the  greater  the  mismatch  in  displacements 
of  adjoining  ceramic  ends.  Since  the  hemispherical 
end  closure  acts  to  support  the  end  of  the  cylinder 
against  buckling,  designing  the  joint  rings  for  addi¬ 
tional  stiffness  would  not  be  considered.  The  first 
assembly  of  figure  3  shows  various  dimensioning 
variables  that  define  the  joint  ring  design.  Selection 
of  appropriate  values  for  these  variables  are  deter¬ 
mined  by  their  impact  on  factors  such  as  hardware 
cost,  ability  to  be  assembled,  and  effect  on  struc¬ 
tural  performance  of  the  joint  interface. 

Selection  of  the  length  of  the  joint  ring  flange  L  is 
driven  by  a  number  of  such  considerations.  The 
longer  the  flange  fi  e.,  the  greater  the  bond  length 
between  cylinder  and  joint  ring),  the  more  expen¬ 
sive  the  part  becomes  to  fabricate.  Other  short¬ 
comings  of  long  flanges  are  the  potential  of 
increased  difficulty  of  assembly  due  to  misalign¬ 
ment  and  a  reduction  in  the  portion  of  the  ceramic 
shell  that  can  be  nondestructively  evaluated  with 
techniques  such  as  ultrasonics.  Intermittent  nonde¬ 
structive  evaluation  of  the  ceramic  housing  wall  to 
ensure  structural  integrity  may  be  desirable  in 
cases  where  the  hull  undergoes  a  high  number  of 
dives  to  design  pressure.  Additionally,  longer- 
length  flanges  may  be  undesirable  from  a  weight 
point  of  view  when  maximum  buoyancy  of  the 
housing  assembly  is  desired. 

The  advantages  of  using  a  longer  flange  are 
primarily  based  on  their  observed  effect  of 
increasing  the  structural  performance  of  the  hous¬ 
ing  assembly  under  cyclic  load.  As  mentioned  pre¬ 
viously,  cracks  that  propagate  from  the  localized 
tensile  stress  regions  &used  by  Poisson’s  effects 
at  the  bearing  surface:  >f  ceramic  ends  can 
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eventually  degrade  the  performance  of  ceramic 
underwater  pressure  housings.  These  circum¬ 
ferential  cracks  have  been  observed  to  run  in  a 
more-or-less  meridional  orientation  from  the  bear¬ 
ing  surface  of  the  ceramic  ends.  Failure  could 
occur  if  these  cracks  propagated  to  the  point  that 
they  break  through  the  inner  or  outer  diameter  of 
the  ceramic  hull  wall  and  a  portion  of  compressive 
load-bearing  hull  spalls  off  or  leaks  occur.  The 
presence  of  a  long  flange  length  L  acts  to  contain 
this  spalling  effect  by  encapsulating  the  ceramic 
ends  where  the  cracks  may  appear.  Additionally, 
finite-element  analysis  (FEA)  indicates  that  the 
presence  of  longer  flange  lengths  also  has  the 
effect  of  decreasing  the  localized  tensile  stresses 
that  may  occur  at  the  ceramic’s  ends  and,  thus, 
reduces  the  chance  of  crack  propagation.  Metallic 
joint  rings  designed  for  ceramic  hemispheres  and 
cylinders  at  NRaD  are  typically  chosen  to  have  a 
flange  length  between  two  and  three  times  the 
thickness  of  the  shell  wall  they  encaps'  'late. 

The  other  primary  variable  that  determines  the 
shape  of  the  flange  is  its  thickness,  TF.  FEA  indi¬ 
cates  that  thinner  flange  thickness  results  in  lower 
tensile  stresses  occurring  at  the  bearing  surfaces 
of  the  ceramic  component  ends.  Manufacturing 
joint  rings  with  thinner  flanges  may  require  tighter 
dimensional  tolerances,  and  thinner  flanges  may 
be  undesirable  if  attachment  points  for  additional 
hardware  are  needed  at  the  flange  surfaces. 

Flange  thickness  may  also  be  constrained  by  the 
amount  of  material  needed  to  contain  the  high 
hydrostatic  stresses  that  occur  in  the  trapped 
epoxy  under  an  external  pressure  load.  The  forces 
exerted  on  the  flanges  of  the  joint  ring  by  the 
epoxy  become  more  critical  the  greater  the  amount 
of  axial  clearance,  TA,  that  exists  between  the 
ceramic  ends  and  the  bearing  surface  on  the  joint 
ring. 

The  amount  of  TA  and  radial  clearance,  TR, 
between  the  ceramic  shell  and  its  metallic  joint 
rings  affects  the  ease  of  joint  ring  assembly  and 
the  structural  performance  of  the  joint  interface. 
The  amount  of  nominal  TR  used  at  NRaD  is  on  the 
order  of  1/1000  of  the  OD  of  the  ceramic  hull  com¬ 
ponents.  Obviously,  the  ability  to  achieve  this 
clearance  in  actual  assembly  depends  on  using 


proper  assembly  fixtures  and  also  on  the  dimen¬ 
sional  tolerances  to  which  the  joint  ring  and 
ceramic  hull  are  fabricated.  Higher  tolerances  on 
the  ceramic  piece  parts  implies  higher  monetary 
costs,  but  also  provides  greater  ease  of  assembly 
and  higher  structural  reliability.  The  additional  cost 
for  tighter  tolerances  is  due  to  the  additional  work 
required  to  finish  grind  the  inner  and  outer  surfaces 
of  the  ceramic  hull.  This  grinding  process  can  be 
used  to  fabricate  ceramic  piece  parts  which  are 
very  concentric  and  of  uniform  wall  thickness.  Typi¬ 
cal  tolerances  for  12-inch  OD  cylinders  (LTD  =  1.5) 
procured  by  NRAD  are  on  the  order  of  plus  or 
minus  .005  of  an  inch  on  wall  thickness  and  OD. 
FEA  indicates  that  using  radial  clearances  of 
1/1 000  of  the  OD  of  the  ceramic  hull  will  result  in 
minimum  interface  stresses  in  the  ceramic  at  the 
joint.  Higher  values  of  radial  clearance  than  this 
could  also  result  in  higher  potential  for  the  epoxy  to 
extrude  from  the  metallic  joint  ring  under  the  nigh 
loads  associated  with  substantial  external  hydro¬ 
static  pressure. 

Controlling  the  TA  between  the  bottom  of  the 
U-shaped  metallic  joint  ring  and  the  bearing  sur¬ 
face  at  the  end  of  the  ceramic  shell  depends  on 
the  relative  flatness  of  these  interfacing  surfaces 
and  the  use  of  standoffs  during  assembly.  A  stand¬ 
off  of  equivalent  thickness  to  the  desired  axial 
clearance  can  be  used  during  assembly  to  main¬ 
tain  spacing  between  the  ceramic  bearing  surface 
and  bottom  of  the  U-shaped  metallic  joint  ring. 
Cured  epoxy  standoffs  or  125-pound  manila  stock 
(see  figure  A-1  of  appendix  A)  can  be  used  for  this 
purpose.  FEA  indicates  that  variations  in  TA  also 
effect  the  interface  stresses  in  the  joint  region.  The 
challenge  of  applying  the  results  of  FEA  for  an 
actual  joint  is  to  develop  manufacturing  tolerances 
and  assembly  techniques  that  achieve  the  desired 
TA  spacing. 

Assembling  the  joint  interface  to  eliminate  epoxy  in 
the  TA  may  also  improve  the  structural  perfor¬ 
mance  of  housings  required  to  achieve  a  high 
number  of  dive  cycles  to  design  depth.  Keeping 
the  bearing  surface  of  the  ceramic  ends  free  of 
epoxy  can  be  accomplished  by  covering  the  bear¬ 
ing  surface  with  a  gasket  (references  6  and  13).  It 
is  hypothesized  that  under  high  hydrostatic  loads 
the  epoxy  could  flow  into  surface  flaws  in  ceramic 
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and  the  high  hydrostatic  pressure  of  the  epoxy  in 
the  flaw  could  then  act  to  initiate  cracks  that  could 
propagate  under  repeated  loading. 

A  number  of  material  candidates  for  use  in  ceramic 
bearing  surface  gaskets  have  been  evaluated 
experimentally  (reference  13)  to  determine  their 
effect  on  the  cyclic  fatigue  life  of  ceramic  underwa¬ 
ter  pressure  housings.  The  most  promising  gasket 
material  tested  based  on  these  studies  is  graphite 
fiber  reinforced  (GFR)  PEEK  composite.  NRaD 
has  used  .040-inch-thick  GFR  PEEK  gaskets  con¬ 
sisting  of  eight  graphite  fiber  plys  laid  up  in  a  (0/90) 
configuration.  In  addition  to  keeping  the  ceramic 
bearing  surface  free  of  epoxy,  it  is  hypothesized 
that  the  tailored  mechanical  properties  of  the  GFR 
PEEK  gasket  are  well  suited  for  application  as  a 
protective  gasket.  The  graphite  fibers  provide  rela¬ 
tively  good  in-plane  stiffness  to  the  gasket  which 
helps  reduce  the  Poisson’s  effect-induced  stresses 
on  the  ceramic  bearing  surface  described  pre¬ 
viously.  Yet,  the  PEEK  matrix  also  provides  the 
gasket  with  tough,  yet  compliant,  properties  in  the 
axial  direction  which  aids  in  reducing  localized 
stresses  on  the  ceramic  bearing  surface  ends  due 
to  surface  irregularities. 

The  benefits  of  having  service  joints  that  can  be 
mechanically  disassembled  are  numerous.  Service 
joints  allow  for  internal  access  to  ceramic  underwa¬ 
ter  pressure  housings  to  package  or  unpackage 
payloads.  They  allow  for  separate  housing  compo¬ 
nents  to  be  pressure  tested  individually  before  the 
entire  assembly  is  tested.  They  allow  for  hull  com¬ 
ponents  to  be  inspected  individually  or  to  be  more 
easily  replaced.  The  use  of  service  joints  also 
requires  additional  sealing  techniques  and  a 
means  of  maintaining  closure.  Figure  3  shows 
three  examples  of  sealing  and  closure  configura¬ 
tions  that  can  be  used  for  service  joints  for  ceramic 
underwater  pressure  housings. 

The  first  assembly  shown  in  figure  3  utilizes  split 
V-band  damp  bands  to  join  adjacent  housing  sec¬ 
tions.  Clamp  bands  aid  in  sealing  a  joint-ring  inter¬ 
face  by  maintaining  uniform  compression  of  O-ring 
face  seals  during  assembly.  These  figures  also 
show  a  redundant  radial  seal  machined  into  the 


extended  lip  of  the  hemisphere  joint  ring.  For  larger 
housings,  pulling  a  slight  internal  vacuum  on  the 
housings  during  assembly  to  compress  the  O-nng 
face  seal  may  help  for  assembling  the  V-band  cou¬ 
pling.  The  location  of  the  O-ring  glands  in  the  joint 
rings  should  be  placed  as  far  as  possible  from  the 
bearing  surfaces  of  the  ceramic  hull  components. 
This  minimizes  the  effect  that  the  potentially  high 
localized  stresses  that  occur  in  the  joint  ring 
around  the  O-ring  gland  under  external  pressure 
loading  will  have  on  the  adjacent  ceramic.  For  this 
reason,  the  O-ring  glands  in  the  joint  rings  shown 
in  figure  3  are  located  in  the  external  flange  of  the 
hemispherical  end-cap  joint  ring.  The  second  and 
third  assemblies  shown  in  figure  3  use  circumfer¬ 
ential  bolts  and  tie  rods  to  maintain  closure.  Ten¬ 
sioned  internal  or  external  tie  rods  between  joint 
rings  of  the  cylindrical  hull  have  the  benefit  of  help¬ 
ing  to  bear  the  handling  loads  on  the  pressure 
housing  assembly  that  would  otherwise  be  bom  by 
the  joint  ring/shell  epoxy  bonds  only.  Using  ten¬ 
sioned  tie  rods  to  precompress  the  ceramic  hull  to 
reduce  tensile  stresses  that  may  occur  in  the 
ceramic  during  handling  is  also  an  option.  The  use 
of  bolts  or  tie  rods  to  maintain  closure  also  requires 
more  space  than  split  V-band  clamp  bands,  which 
may  be  of  concern  if  tight  packaging  volume 
requirements  exist. 

Figure  4  shows  four  detailed  views  of  joint-ring 
designs  based  on  those  discussed  earlier  and 
shown  in  figure  1 .  The  first  two  designs  are  perma¬ 
nent  joints  that  could  be  used  where  disassembly 
is  not  required.  The  first  design  shows  a  coupling 
ring  that  could  be  used  to  assemble  multiple  cylin¬ 
ders  together  to  create  a  longer  hull  for  cases 
when  the  ceramic  cylinder  cannot  be  fabricated  as 
a  single  unit.  The  second  design  shows  a  coupling 
ring  with  an  integral  internal  T-shaped  stiffening 
ring  to  provide  additional  structural  support.  The 
last  two  configurations  shown  in  figure  4  are 
options  when  both  disassembly  and  additional  stiff¬ 
ness  are  required  for  a  joint-ring  design.  The  fourth 
option  has  the  additional  advantage  that  the  inner 
stiffening  l-ring  is  not  externally  exposed,  which 
offers  the  designer  a  greater  selection  of  material 
choices  from  which  to  fabricate  the  1-ring.  Materials 
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which  have  high  specific  modulus,  but  are  suscep¬ 
tible  to  sea  water  corrosion,  can  now  be  consid¬ 
ered. 

Selection  of  the  assembly  tolerances  and  encapsu¬ 
lating  flange  dimensions  discussed  earlier  for  ser¬ 
vice  joint  rings  also  apply  to  the  joints  shown  in 
figure  4.  The  design  of  joint  rings,  like  those  in  the 
final  three  assemblies  of  figure  4,  that  provide  addi¬ 
tional  structural  support  for  buckling  resistance  is 
the  next  issue  to  be  addressed. 

JOINT  RING  DESIGN  FOR  INCREASED 
STRUCTURAL  STABILITY 


Structural  instability  can  occur  in  underwater  pres¬ 
sure  housings  under  external  hydrostatic  load  if 
there  exists  a  means  by  which  the  strain  energy 
associated  with  compression  of  the  shell  mem¬ 
brane  can  be  converted  to  strain  energy 
associated  with  bending  of  the  shell  membrane. 
Typically,  the  in-plane  membrane  stiffness  of  a 
shell  is  substantially  greater  than  its  bending  stiff¬ 
ness.  Consequently,  if  the  strain  energy  in  the 
membrane  can  be  converted  to  bending  energy, 
the  subsequent  bending  of  the  shell  can  result  in 
very  large  deflections.  These  large  deflections 
associated  with  bending  of  the  pressure-housing 
shell  wall  is  known  as  buckling  and  is  one  of  the 
primary  failure  mechanisms  for  ceramic  underwa¬ 
ter  pressure  housings.  Utilizing  stiffening  rings  like 
those  in  the  last  three  assemblies  of  figure  4 
increases  the  buckling  resistance  of  the  housing 
assembly  because  the  inertia  created  by  a  deep 
joint-ring  web  can  result  in  substantial  bending  stiff¬ 
ness. 

The  use  of  a  joint  ring  with  radial  depth  to  provide 
additional  buckling  resistance  is  predicated  on  the 
available  external  and  internal  packaging  volume 
that  exists  for  each  housing.  Optimally,  the  major 
portion  of  mass  of  the  T-shaped  portion  of  a  cen¬ 
tral  stiffened  joint  ring  would  be  concentrated  in  an 
internal  flange  that  is  offset  from  the  shell  wall  via  a 
deep  relatively  slender  web.  External  pressure 
housings  can  also  be  stiffened  with  external  ribs, 
but  this  approach  is  not  as  structurally  efficient  as 
an  internal  rib  configuration.  The  larger  the  diame¬ 


ter  of  the  web  and  flange  of  the  stiffening  ring,  the 
less  stable  it  becomes  under  compressive  load, 
and,  consequently,  the  less  structural  support  it  will 
provide  to  the  hull  assembly. 

Once  volume  constraints  are  established,  stiffened 
joint  rings  require  structural  analysis  techniques  to 
verify  the  integrity  of  the  design.  Stress  analysis 
should  be  performed  to  verify  that  stresses  in  the 
stiffened  joint  ring  are  well  below  levels  that  would 
cause  yielding  of  the  metal  to  occur.  If  yielding 
occurs  before  buckling,  a  stability  analysis  that  can 
account  for  material  nonlinearities  is  required. 
Additionally,  stability  analysis  must  be  performed  to 
ensure  that  the  stiffened  joint  rings  provide  ade¬ 
quate  resistance  against  failure  by  general  instabil¬ 
ity  of  the  housing  (long  wave  length),  as  well  as 
failure  by  local  instabilities  such  as  local  flange  or 
web  crippling  (short  wave  length).  In  designs 
where  multiple  stiffened  joint  rings  are  used,  inner 
bay  buckling  (intermediate  wave  length)  also  must 
be  checked.  An  optimized  hull  design  incorporating 
stiffened  joint  rings  would  have  all  of  these  three 
potential  failure  modes  occur  at  the  same  external 
pressure  load. 

The  amount  of  stiffness  that  a  joint  ring  provides  a 
housing  depends  on  both  its  geometry  and  location 
with  respect  to  the  shell  wall,  as  well  as  the  mate¬ 
rial  that  is  used.  For  this  reason,  high  modulus 
materials  or  high  specific  modulus  materials  where 
weight  is  a  concern  are  most  attractive.  NRaD  has 
utilized  both  titanium  alloy  (Ti-6AI-4V)  and  high- 
strength  7000  series  aluminum  alloys  to  fabricate 
stiffened  joint  rings.  Stiffened  joint  rings  that  utilize 
ceramics  as  part  of  their  design  also  have  been 
considered.  A  stiffened  joint  ring  that  could  utilize 
ceramic's  outstanding  specific  modulus  would  be 
ideal  where  weight  savings  are  of  concern. 

JOINT  RING  TEST  CASES _ 


Four  ceramic  underwater  pressure  housings  were 
assembled  and  tested  to  destruction  to  validate 
potential  joint-ring  designs  for  full-scale  26-  and 
33-inch  OD  ceramic  deep-submergence  pressure 
housings  designed  by  NRaD  (reference  6).  All  four 
test  housings  utilized  alumina-ceramic  cylinders 
purchased  from  COORS  Ceramics  Company  as 
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the  primary  hull  structures.  99.8  percent  alumina 
was  used  for  the  case  1  housing,  while  94  percent 
alumina  was  used  for  the  other  three  housing 
assemblies.  All  metallic  joint  rings  used  in  these 
tests  were  machined  from  Ti-6AI-4V  that  was  sub¬ 
sequently  bonded  to  the  alumina  hull  components 
with  a  two-part  epoxy  as  described  in  appendix  A. 

Joining  two  alumina  cylinders  together  using  a  cen¬ 
tral  coupling  ring  or  a  stiffened  central  joint  ring 
were  two  of  the  concepts  that  were  studied  for  the 
designs  of  the  26-  and  33-inch  housing  cylindrical 
hulls.  The  intent  of  these  four  test  cases  was  to 
compare  their  predicted  performance  based  on 
structural  analysis  calculations  to  their  actual 
tested  performance.  Specifically,  predictions  of  the 
external  pressure  required  for  each  of  the  housings 
to  fail  by  buckling  was  to  be  compared  to  tested 
collapse  pressures,  although  finite-element  stress 
analysis  calculations  were  compared  to  strain  gage 
data  for  test  case  4  (reference  1 3). 

Figure  5  shows  the  final  assembled  housing  used 
for  test  case  1 ,  and  figure  6  shows  an  internal  view 
of  the  case  1  housing  prior  to  bonding  of  the 
second  ceramic  hull  to  the  central  coupling  ring. 
Dimensions  for  each  of  the  ceramic  hulls  used  in 
this  assembly  are  given  in  figure  7,  and  the  dimen¬ 
sions  of  the  titanium  central  coupling  ring  are  pro¬ 
vided  in  figure  8. 

Both  ceramic  shell  components  for  the  case  1 
assembly  have  a  nominal  seven-inch-long  cylindri¬ 
cal  section  with  an  integral  hemisphere  at  one  end. 
A  nominal  wall  thickness  of  .2  of  an  inch  is  used 
for  both  the  cylindrical  and  hemispherical  portions 
of  each  ceramic  hull.  Each  of  the  ceramic  compo¬ 
nents  were  fabricated  by  a  slip  casting  process 
with  no  finishing  operations  performed.  Fabricating 
the  part  without  any  finish  grinding  results  in  a 
lower  cost  part  at  the  expense  of  looser  dimen¬ 
sional  tolera.  ?es.  Loose  dimensional  tolerances 
for  the  ceramic  components  requires  either  loose 
tolerances  for  the  coupling  ring  or  making  the  effort 
to  custom  fit  the  coupling  ring  to  the  measured 
dimensions  of  the  as-fabricated  ceramic  hull  ends. 
Loose  dimensional  tolerances  also  add  more 
uncertainty  to  predicting  failure  by  buckling.  To 


compensate  for  this  increased  uncertainty,  addi¬ 
tional  shell  wall  thickness  may  be  required  for  an 
additional  safety  margin. 

A  cross-sectional  view  of  the  test  case  2  assembly 
is  shown  in  figure  9.  This  configuration  consists  of 
two  alumina-ceramic  cylinders  shown  in  figure  10 
joined  together  using  the  titanium  coupling  ring 
shown  in  figure  11 .  The  remaining  cylinder  ends 
are  encapsulated  with  the  end  cap  joint  rings 
shown  in  figure  1 2  that  are  designed  to  mate  with 
the  flat-steel  end  plates  shown  in  figure  13. 

The  case  3  assembly  is  shown  in  figure  14  and 
utilizes  two  ceramic  cylinders  with  the  same 
dimensions  and  alumina  composition  as  used  for 
the  case  2  assembly  (figure  10).  These  two  cylin¬ 
ders  are  joined  together  using  the  central  stiffened 
joint  ring  shown  in  figure  15  and  their  remaining 
ends  are  encapsulated  with  the  end-cap  joint  rings 
machined  to  the  dimensions  shown  in  figure  12. 
The  central  stiffened  joint  ring  acts  to  couple  the 
two  ceramic  cylinders  together  and  also  provides 
additional  buckling  resistance  through  its  integral 
T-shaped  ring  at  its  inner  wall.  Existing  titanium 
hemispheres  shown  in  figure  1 6  from  an  earlier 
NRaD  program  were  used  as  end  closures  for  the 
case  3  housing.  Figure  1 7  shows  an  end-cap  joint 
ring  and  a  stiffened  central  joint  ring  before  assem¬ 
bly.  Figure  18  shows  each  alumina  cylinder  with  a 
bonded  end-cap  joint  ring  prior  to  epoxy  bonding 
the  stiffened  central  joint  ring  in  place. 

The  configuration  of  the  test  case  4  assembly  is 
covered  in  great  detail  in  reference  13.  The  top 
assembly  drawing  along  with  the  joint  rings  used 
are  shown  again  in  figures  1 9,  20,  21 ,  and  22.  The 
central  stiffened  joint  ring  is  similar  to  the  one  used 
for  test  case  3.  The  short  thick  web  shown  for  the 
integral  T-shaped  ring  does  not  represent  an  opti¬ 
mally  designed  stiffener  because  of  constraints 
that  were  set  for  maximizing  internal  packaging 
volume  while  loading  the  housing  from  one  end 
only.  An  alternative  stiffener  could  be  designed  for 
the  same  weight  as  the  part  shown  in  figure  20 
which  would  provide  substantially  more  external 
pressure  capacity  if  packaging  volume  constraints 
did  not  exist. 
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PROCEDURES  FOR  CALCULATING 
STRUCTURAL  STABILITY  FOR  JOINT 
RING  TEST  CASES 


Structural  stability  analysis  requires  selecting  the 
operating  requirements  (service  depth)  and  desired 
margins  of  safety  for  the  housing  design. 

Traditional  buckling  analysis  involves  calculating 
the  elastic  bifurcation  buckling  pressure  based  on 
perfect  housing  geometry  to  ensure  the  housing 
design  meets  the  selected  requirements.  The  mar¬ 
gin  of  safety  used  for  the  buckling  analysis  is  an 
important  consideration  and  should  depend  largely 
on  the  amount  of  geometric  imperfection  that 
exists  in  the  finished  housing  components.  The 
increased  uncertainty  in  the  buckling  calculations 
brought  about  by  greater  geometric  imperfection 
may  require  more  design  margin  to  ensure  the  per¬ 
formance  on  the  hull  design. 

Preliminary  predictions  of  critical  buckling  pres¬ 
sures  for  cylindrical  pressure  housings  can  be  per¬ 
formed  using  equations  derived  for  closed-ended 
vessels  under  uniform  external  pressure.  The  fol¬ 
lowing  equation  assumes  that  the  ends  of  the 
cylinder  are  simply  supported,  i.e.,  the  ends  are 
rigidly  supported  in  the  radial  direction,  but  meridio¬ 
nal  rotations  are  allowed. 


This  equation  (equation  1), 


(reference  1 ) 

where  n=#  of  lobes  formed,  t=cylinder  wall  thick¬ 
ness,  D=outer  diameter  of  cylinder,  L=simply  sup¬ 
ported  length  of  cylinder,  E=elastic  modulus  of 
cylinder  material,  and  v=Poisson’s  ratio  of  cylinder 
material,  can  be  used  for  initial  predictions  of  col¬ 
lapse  pressure  for  housings  such  as  those  indi¬ 
cated  in  cases  1  and  2.  Derivation  of  equation  1  is 
based  on  linear  differential  equations  where  deflec¬ 
tions  of  the  shell  structure  are  assumed  to  be 
small.  Consequently,  applications  of  this  formula  to 


housings  that  exhibit  nonlinearities  due  to  large 
deformations  and/or  material  nonlinearities  would 
result  in  less  accurate  predictions  of  failure  by 
buckling. 

A  solid  flat-end  plate  with  a  short  lap  interface  with 
the  cylinder,  such  as  used  in  case  2,  is  reasonably 
approximated  to  be  a  simply  supported  boundary 
condition.  In  cases  where  the  cylinder  is  capped  at 
both  ends  by  hemispheres,  equation  1  can  still  be 
employed  using  an  effective  length  (L)  equal  to  the 
simply  supported  length  of  the  cylinder,  plus  one 
third  the  depth  of  each  hemisphere.  Trie  logic  *or 
using  an  effective  length  follows  since  a  hemi¬ 
sphere  represents  a  more  compliant  end  constraint 
than  a  flat  plate  and,  thus,  would  be  expected  to 
decrease  the  collapse  pressure  for  a  cylindrical 
housing  in  the  same  way  that  lengthening  the 
cylindrical  portion  of  a  housing  would  also 
decrease  the  collapse  pressure.  Further,  it  has 
been  shown  (reference  8)  that  varying  the  wall 
thickness  of  the  hemisphere  or  using  a  different 
modulus  material  for  the  hemisphere  does  rela¬ 
tively  little  to  change  the  collapse  pressure  of  trie 
housing  assembly.  Thus,  increasing  the  effective 
length  (L)  of  the  housing  by  one  third  the  depth  of 
each  hemisphere  regardless  of  hemisphere  design 
details  is  a  reasonable  first  approach. 

Since  the  hemispheres  in  case  1  are  integral  end 
closures  to  the  cylinders,  the  use  of  equation  1 
which  assumes  simply  supported  boundary  condi¬ 
tions  at  cylinder  ends  would  be  expected  to  give 
conservative  results.  Yet,  the  presence  of  addi¬ 
tional  boundary  conditions  at  the  cylinder  ends  has 
relatively  little  influence  on  collapse  pressures  due 
to  general  instability  failure  modes  for  the  entire 
cylindrical  housing  assembly  unless  the  length  of 
the  cylindrical  section  is  relatively  short. 

Equation  1  also  can  be  used  when  separate  cylin¬ 
drical  sections  are  joined  with  coupling  rings  such 
as  in  cases  1  and  2  by  assuming  that  the  cylindri¬ 
cal  portion  of  the  housing  behaves  as  a  monoco- 
que  cylinder  with  length  (L)  equal  to  sum  of  the 
length  of  the  two  cylinders  and  the  web  thickness 
of  the  coupling  ring.  For  the  coupling  rings  used  in 
cases  1  and  2,  this  has  been  shown  to  be  a  con¬ 
servative  means  of  calculating  the  collapse  pres¬ 
sure  since  the  increased  inertia  associated  with  the 
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flanges  of  the  coupling  ring  actually  help  to 
increase  buckling  resistance  of  the  housing.  While 
equation  1  may  be  used  to  help  predict  the  general 
instability  of  housings  joined  by  coupling  rings, 
other  techniques  need  to  be  used  to  ensure  that 
the  flanges  of  the  coupling  ring  are  substantial 
enough  to  resist  local  crippling  under  external  pres¬ 
sure. 

The  increase  in  buckling  resistance  afforded  by 
adding  more  material  to  the  flanges  of  the  coupling 
ring  motivates  the  design  of  the  joint  rings  used  in 
cases  3  and  4.  If  the  stiffness  of  the  joint  ring  was 
substantial  enough  to  approximate  its  support  of 
the  cylinder  as  being  simply  supported,  then  equa¬ 
tion  1  could  be  used  to  predict  innerbay  buckling 
for  each  of  the  cylinders.  For  the  stiffened  central 
joint  rings  used  in  cases  3  and  4,  the  critical  col¬ 
lapse  mode  is  by  general  instability  for  the  entire 
cylindrical  assembly  (two  cylinders  joined  by  a  stiff¬ 
ened  central  joint  ring).  Hand  calculation  tech¬ 
niques  for  predicting  failure  for  this  type  of 
assembly  are  limited  such  that  stability  analysis 
requires  the  use  of  more  comprehensive  tools 
such  as  the  BOSOR4  computer  program. 

BOSOR4  is  a  structural  analysis  program  for  com¬ 
puting  stress,  buckling,  and  vibration  of  complex 
shells  of  revolution  that  was  developed  by  David 
Bushnell  at  Lockheed  Missiles  and  Space  Co.,  Inc. 
BOSOR4  is  written  using  FORTRAN  IV,  and  pre¬ 
dictions  are  based  on  finite  difference  energy  mini¬ 
mization  with  constraint  conditions  for  the  structure 
under  study.  The  meridian  of  the  shell  of  revolution 
is  modeled  by  using  a  number  of  segments  with 
material,  geometric,  and  boundary  condition  prop¬ 
erties  representative  of  the  real  structure. 

B0S0R4  predictions  of  collapse  pressure  for  the 
four  test  cases  presented  here  were  based  on 
elastic  bifurcation  analysis.  All  housing  materials 
were  modeled  as  linear  elastic,  but  nonlinear 
effects  due  to  large  deflections  were  considered. 
Bifurcation  analysis  searches  for  the  external  pres¬ 
sure  load  at  which  the  axisymmetric  deformation  of 
the  housing  structure  ceases  to  be  stable.  For 
externally  pressurized  cylindrical  housings  like  the 


four  test-case  assemblies  presented  here,  the  post 
buckling  behavior  is  unstable  such  that  once  the 
bifurcation  pressure  is  reached,  the  housing  can¬ 
not  continue  to  carry  additional  external  pressure 
which  results  in  catastrophic  collapse  of  the  hull. 

For  real  housing  structures,  there  is  no  such  thing 
as  true  bifurcation  buckling.  All  real  housings  have 
imperfections  which  lead  to  discrepencies  between 
predicted  failure  by  bifurcation  theory  and  the 
actual  tested  failure.  The  amount  of  discrepency  is 
dependent  on  the  amount  of  initial  imperfection  in 
the  housing  components.  Cylindrical  housings  sub¬ 
jected  to  external  pressure  are  imperfection  sensi¬ 
tive,  meaning  the  presence  of  initial  imperfections 
act  to  reduce  the  pressure  capacity  of  the  housing 
as  predicted  by  bifurcation  theory.  Cylindrical  hous¬ 
ings  subjected  to  hydrostatic  external  pressure 
typically  have  general  instability  failure  modes 
characterized  by  long  axial  wavelengths  and  rela¬ 
tively  few  circumferential  wavelengths  (N=2  or  N=3 
for  test  cases  discussed  in  this  report).  Conse¬ 
quently,  global  imperfections  such  as  uniform  out- 
of-roundness  would  be  expected  to  have  more  del¬ 
eterious  effects  than  any  localized  imperfections. 

A  linear  elastic  bifurcation  analysis  of  externally 
pressurized  housings  is  appropriate  when  there  is 
relatively  little  bending  energy  stored  in  the  hous¬ 
ing  up  until  the  bifurcation  point  is  reached.  This  is 
the  case  in  housings  where  there  is  a  good  radial 
deflection  match  between  adjacent  housing  com¬ 
ponents  implying  the  axisym metrically  deformed 
housing  has  the  same  approximate  shape  as  the 
undeformed  housing  prior  to  loading.  Conversely,  if 
substantial  bending  energy  exists  in  the  cylindrical 
housing  subjected  to  external  pressure  due  to  rela¬ 
tively  rigid  boundary  conditions  such  as  occurs  at 
the  interface  with  flat  bulkhead  end  closures  or  at 
stiff  joint  rings,  then  a  nonlinear  analysis  may  be 
required.  The  power  of  programs  like  BOSOR4  is 
that  they  allow  the  designer  of  underwater  pres¬ 
sure  housings  to  perform  buckling  analysis  based 
on  different  modeling  assumptions  (i.e.,  linear  or 
nonlinear  geometric  or  material  effects)  as 
appropriate  for  each  particular  housing  design 
case. 
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RESULTS  OF  JOINT  RING 
STRUCTURAL  STABILITY 
CALCULATIONS 


The  buckling  capacity  of  the  four  test-case  housing 
assemblies  were  predicted  using  BOSOR4  and 
hand  calculation  techniques  where  applicable.  The 
results  of  these  calculations  are  shown  in  table  1 . 
The  material  properties  for  the  ceramic  cylinders 
used  for  each  housing  assembly  are  listed  in  the 
table  along  with  the  housing  dimensions  used  for 
calculating  collapse  pressure  for  cases  1  and  2 
using  equation  1  shown  above.  All  titanium  joint 
rings  were  modeled  using  material  property  values 
of  16.4  million  for  the  compressive  elastic  modulus 
and  .31  for  the  Poisson’s  ratio.  Equation  1  was 
programmed  in  FORTRAN  to  calculate  collapse 
pressure  of  the  first  two  test  cases  for  values  of 
N=2,  3,  and  4.  The  results  of  pressure  tests  per¬ 
formed  with  each  of  the  housings  are  listed  for 
comparison  at  the  bottom  of  table  1 . 

The  BOSOR4  models  were  constructed  to  allow 
for  meridional  rotations  at  the  joints  between  cylin¬ 
ders  and  end  closures  under  external  pressure  as 
occurs  in  cases  2  through  4.  The  central  coupling 
rings  and  central  stiffened  joint  ring  were  modeled 
with  their  own  discrete  branched-shell  segments. 
Detailed  modeling  of  these  central  joint  rings 
allows  for  prediction  of  any  local  instabilities  in  the 
ring  as  well  as  more  accurate  predictions  of  gen¬ 
eral  instability  modes  for  the  entire  housing. 

Results  of  hand  calculations  show  prediction  of  the 
collapse  pressure  based  on  equation  1  to  be 
11,582  psi,  which  is  off  by  more  than  20  percent 
from  the  actual  tested  value  of  9,200  psi  for 
case  1 .  The  simply  supported  length  (L)  used  in 
equation  1  for  the  case  1  housing  was  16.16 
inches  (7.00  inches  for  each  ceramic  cylinder  sec¬ 
tion,  plus  1 .02  inches  for  each  ceramic  hemisphere 
section  equal  to  one  third  the  depth  of  hemisphere, 
and  .12  of  an  inch  for  coupling  ring  web).  The 
BOSOR4  model  for  the  case  1  housing  consisted 
of  seven  segments.  Two  segments  were  used  to 
model  each  ceramic  hull  component,  and  a  total  of 
three  segments  were  used  to  model  the  web  and 
two  flanges  of  the  central  coupling  ring.  BOSOR4 
calculations  for  the  case  1  housing  predicted  a 


collapse  pressure  of  10,512  psi,  which  also  was 
higher  than  the  actual  tested  value.  Some  of  this 
difference  could  be  attributed  to  the  deleterious 
effect  of  geometric  imperfections  that  existed  in  the 
ceramic  hull  components  that  were  used. 

The  ceramic  components  used  in  the  test  case  1 
housing  had  maximum  deviations  from  true  circular 
shape  equal  to  about  10  percent  of  the  nominal 
shell  wall  thickness.  BOSOR4  models  for  all  four 
test  cases  were  run  using  nominal  housing  compo¬ 
nent  dimensions  and  the  assumption  that  no  geo¬ 
metric  imperfections  were  present.  Additionally,  the 
use  of  equation  1  assumes  meridional  rotations 
can  occur  between  the  end  closure  and  the  cylin¬ 
drical  hull  which  does  not  hold  for  case  1  since  the 
hemispheres  are  integral  to  the  cylinder  sections  of 
the  housing. 

Both  hand  calculations  and  the  BOSOR4  model 
predict  that  the  case  1  housing  would  fail  with  two 
circumferential  lobes  forming  as  the  housing 
buckled.  Figure  23  shows  BOSOR4  graphical  out¬ 
put  for  the  predicted  failure  mode  for  case  1 .  This 
figure  shows  the  meridian  of  the  shell  before  pres¬ 
sure  loading  (dashed  lines)  and  its  buckled  config¬ 
uration  as  it  fails  (solid  lines).  Obviously,  the 
formation  of  lobes  is  not  an  axisymmetric  deforma¬ 
tion  such  that  the  buckled  configuration  shown  in 
figure  23  represents  a  sectioned  view  through  the 
major  axis  of  a  single  circumferential  wave.  As  a 
comparison,  a  BOSOR4  model  was  constructed 
that  contained  a  single  monocoque  cylinder  with  an 
equivalent  length  to  the  two-cylinder/coupling-ring 
configuration  of  case  1  (outer  diameter  and  wall 
thickness  of  the  ceramic  monocoque  cylinder  was 
kept  constant).  Collapse  by  buckling  was  predicted 
to  occur  at  10.201  psi  for  this  single-cylinder 
replacement  which  indicates  the  presence  of  this 
particular  coupling  ring  slightly  increases  the  exter¬ 
nal  pressure  capacity  of  the  case  1  assembly 
(10,201  psi  to  10,512  psi).  This  indicates  that  a 
cylindrical  hull  could  be  built  up  using  a  number  of 
cylinders  joined  by  coupling  rings  without  neces¬ 
sarily  degrading  the  structural  performance  of  the 
cylindrical  hull  assembly. 

Calculations  based  on  equation  1  and  a  BOSOR4 
model  both  under-predict  the  buckling  capacity  for 
the  case  2  housing,  although  both  values  are  still 
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within  a  reasonable  margin  of  the  tested  value. 

Use  of  equation  1  predicts  a  collapse  pressure  of 
8,413  psi,  and  a  BOSOR4  model  predicts  a  col¬ 
lapse  pressure  of  9,558  psi  as  compared  to  the 
tested  result  of  a  critical  pressure  of  10,250  psi. 
The  simply  supported  length  (L)  used  in  equation  1 
for  the  case  2  housing  was  29.03  inches 
(1 5.00  inches  for  each  cylinder,  plus  .25  of  an  inch 
for  coupling  ring  web,  and  minus  .61  of  an  inch  for 
supported  length  at  each  flat-steel  end  plate).  The 
model  used  for  the  case  2  housing  consisted  of 
seven  segments.  One  segment  was  used  to  model 
each  end  plate,  and  one  segment  was  used  to 
model  each  ceramic  cylinder.  Three  segments 
were  used  to  model  the  web  and  two  flanges  of  the 
central  coupling  ring.  Initially,  additional  segments 
were  used  to  model  both  end-cap  joint  rings,  but 
they  were  eliminated  in  the  final  model  since  their 
presence  had  virtually  no  effect  on  the  general 
instability  modes  for  the  entire  assembly.  The 
BOSOR4  model  for  the  case  2  housing  accounts 
for  the  increased  support  associated  with  the  stiff¬ 
ness  supplied  by  the  flanges  of  the  coupling  ring 
and  gives  a  more  accurate  prediction  than 
obtained  with  equation  1  where  the  coupling  ring 
flanges  are  unaccounted  for. 

Both  techniques  predict  that  failure  should  occur 
for  the  case  2  assembly  with  three  circumferential 
waves  occurring  at  collapse.  Figure  24  shows  the 
BOSOR4  graphical  output  for  the  predicted  failure 
mode  for  the  case  2  assembly.  Figure  25  shows 
the  case  2  coupling  ring  after  the  test  housing  was 
pressurized  to  failure.  Visual  inspection  of  the 
deformed  coupling  ring  indicated  that  it  was  slightly 
triangular  in  shape  as  would  be  expected  for  an 
N=3  failure.  The  eventual  formation  of  circumferen¬ 
tial  lobes  results  in  bending  the  housing  shell  wall 
to  the  point  that  rupture  occurs.  In  figure  25,  the 
external  flange  of  the  coupling  ring  can  be  seen  to 
have  blown  outward  at  an  assumed  location  of  one 
of  the  three  circumferential  lobes  that  formed  as 
the  housing  failed.  Figure  26  shows  the  end-cap 
joint  rings  from  the  failed  case  2  test  housing. 
These  joint  rings  showed  substantially  less  perma¬ 
nent  deformation  than  occurred  with  the  central 
coupling  ring  as  would  be  expected  since  the  pre¬ 
dicted  maximum  deflections  for  the  case  2 


assembly  should  occur  at  the  mid-bay  location  of 
the  central  coupling  ring  prior  to  failure. 

Buckling  predictions  were  not  made  with  equa¬ 
tion  1  for  cases  3  and  4  because  of  the  complicat¬ 
ing  presence  of  the  central  stiffened  joint  rings 
utilized  in  both  these  assemblies.  However,  predic¬ 
tions  were  made  based  on  BOSOR4  models  for 
these  case  3  and  case  4  housings  which  gave 
excellent  results.  The  case  3  housing  was  pre¬ 
dicted  to  fail  with  two  lobes  forming  (N=2)  at  an 
external  pressure  of  9,603  psi  compared  with  the 
tested  failure  pressure  of  9,550  psi.  The  BOSOR4 
model  for  this  test  case  used  a  toted  of  nine  seg¬ 
ments.  One  segment  was  used  to  model  each  tita¬ 
nium  hemisphere,  and  one  segment  was  used  to 
model  each  ceramic  cylinder.  A  total  of  five  seg¬ 
ments  was  used  to  model  the  three  flanges  and 
two  webs  that  make  up  the  central  stiffened  joint 
ring. 

Figure  27  shows  BOSOR4  graphical  output  for  the 
predicted  mode  of  failure  for  case  3.  Figure  28 
shows  the  hardware  that  was  retrieved  after  pres¬ 
sure  testing  the  housing  to  failure.  The  remains  of 
a  wood  plug  used  to  mitigate  the  effects  of  the 
implosion  on  the  titanium  hemispheres  can  be 
seen.  Figure  29  shows  the  remnants  of  the  central 
stiffened  coupling  ring  which  clearly  validates  the 
prediction  that  failure  occurred  with  two  circumfer¬ 
ential  lobes  forming  due  to  an  N=2  general  instabil¬ 
ity.  It  is  apparent  from  this  picture  that  the  wooden 
implosion  mitigation  plug  did  relatively  little  to  resist 
snap-through  of  the  stiffened  central  joint  ring  as 
the  housing  assembly  buckled. 

Detailed  BOSOR4  modeling  of  the  geometry  of  the 
central  stiffened  joint  ring  also  allows  the  designer 
to  predict  any  buckling  failures  that  may  occur  due 
to  local  instabilities  such  as  crippling  of  joint  ring 
flanges  or  webs.  BOSOR4  calculations  predicted 
that  local  buckling  of  the  central  stiffened  joint  ring 
for  case  3  would  occur  at  external  pressures  sub¬ 
stantially  greater  than  that  required  to  initiate  fail¬ 
ure  by  general  instability.  Ideally,  the  design  of  the 
central  stiffened  joint  ring  would  be  optimized  such 
t^at  local  crippling  would  occur  at  the  same  exter¬ 
nal  pressure  that  causes  a  general  instability 
failure.  But  as  mentioned  earlier,  packaging 
constraints  limited  the  geometry  of  the  central 
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stiffened  joint  ring  and,  thus,  prevented  an  opti¬ 
mized  design. 

Prior  to  pressurizing  the  case  3  housing  to  failure, 
the  housing  assembly  was  equipped  with  strain 
gages  on  the  interior  and  exterior  surfaces  at  mid¬ 
bay  of  each  ceramic  cylinder  as  well  as  on  various 
locations  of  the  central  stiffened  joint  ring  as  indi¬ 
cated  in  figure  30.  Varying  numbers  of  gages  were 
located  around  the  circumference  of  the  housing 
for  each  axial  location  shown.  Eleven  pressure 
cycles  to  9,000  psi  were  completed,  and  strains 
were  recorded  for  each  gage  at  1 ,000  psi  intervals 
until  an  external  pressure  of  8,500  psi  was  reached 
whereby  the  recording  interval  was  reduced  to 
1 00  psi  increments.  Table  2,  Sheets  1  through  6 
show  strain  gage  measurements  recorded  during 
testing  in  units  of  micro  inches/inch.  Plots  of  th.s 
data  are  presented  in  figures  31  through  34.  These 
plots  reveal  that  the  measured  strains  increase  as 
a  linear  function  of  pressure  until  an  external  pres¬ 
sure  of  8,500  psi  is  reached.  At  external  pressures 
greater  than  8,500  psi,  the  strains  recorded  for  the 
various  circumferential  gages  at  each  axial  location 
begin  to  diverge. 

The  initiation  of  this  nonlinear  behavior  indicates 
that  circumferential  waves  are  forming.  It  is  of 
value  to  note  that  the  onset  of  lobes  occurs  well 
before  the  housing  assembly  finally  collapses. 
Specifically  for  the  case  3  housing,  lobes  have 
begun  to  form  at  an  external  pressure  equal  to 
1 ,000  psi  less  than  the  collapse  pressure.  As 
expected,  the  formation  of  circumferential  waves  is 
most  accentuated  in  the  data  for  the  inner  flange  of 
the  central  stiffened  joint  ring  shown  in  figure  33. 
The  additional  bending  of  the  shell  associated  with 
the  formation  of  lobes  eithe-  acts  to  increase  or 
reduce  the  measured  strains  depending  upon  the 
location  of  the  gage  with  respect  to  the  waves  that 
form.  It  is  this  bending  which  accounts  for  the  non- 
linearities  that  appear  in  the  data  near  the  collapse 
pressure  for  this  particular  housing  assembly. 

The  BOSOR4  prediction  of  a  critical  buckling  pres¬ 
sure  of  12,208  psi  for  the  case  4  housing  assembly 
compared  very  well  with  the  tested  collapse  pres¬ 
sure  of  11 ,930  psi.  The  BOSOR4  model  for  the 
case  4  housing  utilized  a  total  of  1 5  segments  to 
define  the  complete  assembly.  Three  segments 


were  used  to  model  each  hemisphere  to  account 
for  the  spherical,  conical,  and  cylindrical  portions  of 
the  hemisphere  design.  One  segment  was  used  to 
model  each  ceramic  cylinder,  and  seven  segments 
were  used  to  model  the  central  stiffened  joint  ring. 
Again,  additional  segments  were  used  to  model 
the  cylinder  and  hemisphere  end-cap  joint  rings  to 
check  for  local  instabilities,  but  were  removed  for 
the  final  model  because  of  their  negligible  impact 
on  the  general  instability  failure  modes  for  the 
whole  assembly.  Appendix  B  contains  the 
BOSOR4  input  file  for  the  case  4  housing  for  refer¬ 
ence  of  the  modeling  techniques  used  for  all  four 
test  cases  presented  in  this  report. 

Figure  35  shows  the  BOSOR4  graphical  output  for 
the  N=2  predicted  failure  for  this  particular  hous¬ 
ing.  Figure  36  shows  the  case  4  housing  as  it  is 
being  loaded  into  the  vessel  for  pressurization  to 
failure.  Figure  37  shows  the  central  stiffened  joint 
ring  after  collapse.  The  case  4  housing  was  par¬ 
tially  filled  with  water  prior  to  pressurizing  to  failure 
to  mitigate  the  effects  of  the  implosion.  This 
resulted  in  substantially  less  deformation  in  the 
central  stiffened  joint  ring  than  in  the  case  3  hous¬ 
ing.  Subsequent  measurements  determined  that 
the  joint  ring  was  oval  in  shape,  corresponding  to 
the  formation  of  two  lobes  (N=2  failure). 

The  N=3  failure  mode  for  the  case  3  and  4  hous¬ 
ings  occur  at  external  pressures  substantially 
higher  than  the  pressure  required  for  the  housings 
to  fail  in  an  N=2  mode.  None  the  less,  it  is  interest¬ 
ing  to  note  the  differences  in  the  way  the  housings 
fail  for  these  two  modes.  The  N=3  failure  mode 
occurs  with  antisymmetric  inner-bay  lobes  forming 
where  the  three  lobes  formed  in  one  cylinder  are 
out  of  phase  with  the  three  lobes  formed  in  the 
adjacent  cylinder  as  shown  in  figure  38.  For  this 
N=3  mode,  the  stiff  central  coupling  ring  merely 
twists  and  has  little  impact  on  the  collapse  pres¬ 
sure,  and  the  shell  wall  thickness  of  the  cylinders 
controls  the  buckling  capacity.  This  type  of  failure 
mode  is  referred  to  as  “rolling”  and  is  caused  by 
compression  of  the  web  of  the  internally  central 
stiffened  joint  ring.  This  implies  that  an  optimal 
housing  consisting  of  two  cylinders  and  a  stiff  cen¬ 
tral  coupling  ring  can  be  designed  such  that  failure 
pressure  would  be  equal  for  both  the  N=2  and  N=3 
mode  shapes.  The  central  coupling  ring  would  be 
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sized  to  provide  adequate  bending  resistance 
against  N=2  failure,  and  the  cylinder  wall  thickness 
would  be  sized  to  be  adequate  for  N=3  failure.  Of 
course,  the  final  design  of  the  coupling  ring  and 
cylinders  also  must  be  satisfactory  from  a  stress 
analysis  point  of  view  and  checked  for  any  local 
instabilities. 

As  an  interesting  comparison,  a  BOSOR4  model 
was  constructed  where  the  case  4  cylindrical  hull 
consisting  of  two  ceramic  cylinders  with  stiffened 
central  joint  rings  was  replaced  by  a  single  mono- 
coque  cylinder  of  equivalent  length  (cylinder  wall 
thickness  and  outer  diameter  were  held  constant). 
This  new  cylindrical  hull  with  case  4  end  closures 
was  predicted  to  buckle  at  8,234  psi.  Thus,  the 
integration  of  the  case  4  central  stiffened  joint  ring 
into  the  cylindrical  hull  assembly  increases  the 
buckling  capacity  by  approximately  48  percent 
(8,234  psi  to  12,208  psi). 


CONCLUSION 


The  use  of  epoxy-bonded  metallic  joint  rings  to 
assemble  ceramic  hull  sections  has  proven  to  be 
an  economically  viable  as  well  as  structurally  reli¬ 
able  means  of  constructing  ceramic  underwater 
pressure  housings  without  size  limitations.  Accu¬ 
rate  predictions  of  the  structural  performance  of 
composite  housings  consisting  of  ceramic  shells 
joined  by  metallic  joint  rings  can  be  made  using 
both  commercially  available  structural  analysis 
software  as  well  as  hand  calculation  techniques. 
Metallic  joint  rings  also  provide  engineers  with  a 
way  to  integrate  design  features  such  as  local 
attachment  points,  additional  buckling  resistance, 
and  a  means  of  sealing  and  maintaining  closure 
between  adjacent  housing  sections.  The  continued 
development  of  simple  joining  techniques  like 
epoxy-bonded  metallic  joint  rings  is  essential  to  the 
continued  success  and  acceptance  of  ceramic 
underwater  pressure  housings. 
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GLOSSARY 


AUTB 

Advanced  Ultrasonic  Test  Bed 

FEA 

FEM 

finite  element  analysis 
finite  element  model 

GFR 

GFRP 

graphite-fiber-reinforced 
graphite-fiber-reinforced  plastic 

ID 

IED 

inner  diameter 

independent  exploratory  development 

Kpsi 

one  thousand  psi 

L 

L/D 

length 

length/diameter 

MEK 

MOR 

methyl  ethyl  ketone 

Modulus  of  rupture 

ND 

nondestructive 

NDE 

nondestructive  evaluation 

NDT 

nondestructive  test 

NOSC 

Naval  Ocean  Systems  Center 

OD 

outside  diameter 

PEEK 

poly-ether-ether-ketone 

rms 

root  mean  square 

SAM 

Scanning  Acoustic  Microscopy 

S.F. 

safety  factor 

t 

ceramic  shell  thickness 

TA 

axial  clearance 

t/D 

thickness/diameter 

TF 

flange  thickness 

TR 

radial  clearance 

t/Ro 

thickness/extemal  radius 

W 

width 

W/D 

weight-to-displacement 

17 


FEATURED  RESEARCH 


18 


FEATURED  RESEARCH 


REFERENCES 


1.  Bickell,  M.  B.  and  Dr.  M.  Ruiz.  “Pressure  Ves¬ 
sel  Design  and  Analysis,”  St.  Martin’s  Press. 

2.  Brush,  Don  O.  and  Bo  O.  Almroth.  1975. 
“Buckling  of  Bars,  Plates,  and  Shells,” 
McGraw-Hill  Company. 

3.  Bushnell,  David.  “BOSOR4  User’s  Manual,” 
Lockheed  Missiles  &  Space  Co.,  Inc. 

4.  Bushnell,  David.  1981.  “Buckling  of  Shells-Pit- 
fall  for  Designers,”  AIAA  Journal,  Vol.  19 
(Sep),  pg.  1183. 

5.  Bushnell,  David.  1981.  “Computerized  Buck¬ 
ling  Analysis  of  Shells,”  Applied  Mechanics 
Laboratory,  Lockheed  Palo  Alto  Research 
Laboratory  (Dec). 

6.  Johnson,  R.  P.,  R.  R.  Kurkchubascne,  and  J. 
D.  Stachiw.  1993.  “Design  and  Structural 
Analysis  of  Alumina  Ceramic  Housings  for 
Deep  Submergence  Service  Fifth  Generation 
Housings,”  NRaD  TR  1583  (Mar),  NCCOSC 
RDT&E  Division,  San  Diego,  CA. 

7.  Johnson,  R.  P.  “Stress  Analysis  Consider¬ 
ations  for  Deep  Submergence  Ceramic  Pres¬ 
sure  Housings,”  ROV  92. 

8.  Kurkchubasche,  R.  R.  “Elastic  Stability  Con- 
sideratins  for  Deep  Submergence  Ceramic 
Pressure  Housings,”  ROV  92. 


9.  Santella,  M.  L  1992.  “A  Review  of  Techniques 
for  Joint  Advance  Ceramics,"  American 
Ceramic  Society  Bulletin,  Vol.  71,  No.  6. 

10.  Stachiw,  J.  D.  1990.  “Exploratory  Evaluation 
of  Alumina  Ceramic  Housings  for  Deep  Sub¬ 
mergence  Service — Fourth  Generation  Hous¬ 
ings,”  NOSC  TR  1355  (Jun),  Naval  Ocean 
Systems  Center,  San  Diego,  CA. 

11.  Stachiw,  J.  D.,  and  J.  L  Held.  1987  “Explor¬ 
atory  Evaluation  of  Alumina  Veramic  Cylindri¬ 
cal  Housings  for  Deep  Submergence 
Service — The  Second  Generation  NOSC 
Ceramic  Housings,”  NOSC  TR  1176  (Sep), 
Naval  Ocean  Systems  Center,  San  Diego, 

CA. 

12.  Stachiw,  J.  D.,  T.  J.  Henderson,  and  C.  A. 
Andersson.  1991.  “Novel  Ceramic  Matrix 
Composites  for  Deep  Submergence  Pressure 
Vessel  Applications,"  NOSC  TD  2222  (Oct), 
Naval  Ocean  Systems  Center,  San  Diego, 

CA. 

13.  Stachiw,  J.  D„  R.  P.  Johnson,  and  R.  R.  Kurk¬ 
chubasche.  1993.  “Evaluation  of  Scale  Model 
Ceramic  Pressure  Housing  for  Deep  Submer¬ 
gence  Service,”  NRaD  TR  1582  (Mar), 
NCCOSC  RDT&E  Division,  San  Diego,  CA. 

14.  Young,  Warren  C.  1989.  “Roark’s  Formulas 
for  Stress  and  Strains,  Sixth  Edition,” 
McGraw-Hill  Company. 


19 


FEATURED  RESEARCH 


METALLIC  END  CAP 


CERAMIC  CYLINDRICAL  HULL 
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Figure  1 .  Three  potential  cylindrical  hull  assemblies  for  cer-vm  underwater  pressure 
housings. 
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Figure  2.  Epoxy  bond  between  joint  ring  and  ceramic  hull. 
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STIFFENED  JOINT  RING 
MECHANICALLY  ASSEMBLED 


Figure  4.  Stiffened  joint  ring  designs. 
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Figure  7.  Case  1  ceramic  hull  component  details. 
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FLAT  STEEL  END  CLOSURE 


CASE  2  COUPLING  RING 
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SK9402-111 
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flat  steel  end  closure  details. 
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Figure  1 5.  Case  3  central  stiffened  joint  ring  details. 
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Figure  16.  Case  3  hemispherical  end  closure  details. 
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Figure  19.  Case  4  test  assembly  for  evaluation  of  central  stiffened  joint  rings,  Sheet  1 
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9.  Case  4  test  assembly  for  evaluation  of  central  stiffened  joint  rings,  Sheet  2. 
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Figure  20.  Case  4  central  stiffened  joint  ring  details. 
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hemisphere  end  cap  joint  ring  details. 


Figure  22.  Case  4  cylinder  end  cap  joint  ring  details. 
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Figure  24.  BOSOR4  predicted  failure  mode  for  case  2  test  assembly. 
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Figure  28.  Case  3  test  assembly  after  failure  by  buckling. 
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LOCATION  OF  STRAIN  GAGES 
ON 

CENTRAL  JOINT  STIFFENER 


STRAIN  GAGES:  CEA-06-125WT-120 

Figure  30.  Location  of  strain  gages  on  case  3  central  stiffened  joint  ring. 
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Figure  31 .  Plot  of  strains  recorded  on  interior  surface  of  case  3  ceramic  cylinders  at  midbay. 
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Figure  33.  Plot  of  strains  recorded  on  interior  surfaces  of  case  3  central  stiffened  joint  ring. 
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Figure  34.  Plot  of  strains  recorded  on  exterior  surfaces  of  case  3  central  stiffened  joint  ring. 
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Figure  35.  B0S0R4  predicted  failure  mode  for  case  4  test  assembly. 
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Figure  37.  Buckle .1  case  4  central  stiffened  joint  ring. 
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Table  1 .  Results  of  joint  ring  structural  stability  calculations. 


CASE  1 

CASE  2 

CASE  3 

CASE  4 

E=54E6  psi 
v=.22 
L=16.16" 

D=6.1’ 

t=.2* 

E=44E6  psi 
v=,21 
L=29.03” 
D=11.89” 
t=.355* 

E=44E6  psi 
v=.21 

E=44E6  psi 
v=.21 

EQUATION  1 

Pent  (psi) 

Pcrit  (psi) 

N=2 

11582 

9274 

N=3 

13253 

8413 

N=4 

22141 

13680 

B0S0R4 

Pcrit  (psi) 

Pcrit  (psi) 

Pcrit  (psi) 

Pcrit  (psi) 

N=2 

10512 

13083 

12208 

N=3 

13978 

9558 

11771 

18687 

N=4 

24507 

15245 

15338 

26919 

TEST 

RESULTS 

BUCKLED  @ 
9200  psi 

BUCKLED  @ 
10,250  psi 

BUCKLED  @ 

9550  psi 

BUCKLED  @ 
11,930  psi 
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Table  2.  Strains  recorded  during  pressure  testing  of  case  3  test  assembly,  Sheet  1 . 


|< - - Cylinder  Ho.  1  Midbay  Interior . . -»| 


Pressure 

Hoop 

Axial 

Hoop 

Hoop 

Hoop 

Hoop 

Axial 

Hoop 

Hoop 

Hoop 

Hoop 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1000 

-350 

-103 

-351 

-354 

-3  .' 

-349 

-106 

-353 

-354 

-355 

-354 

2000 

-692 

-213 

-695 

-699 

-696 

-692 

-219 

-696 

•699 

-698 

-694 

3000 

-1040 

-326 

-1046 

-1051 

-1045 

-1041 

-330 

-1046 

-1051 

-1049 

-1040 

4000 

-1380 

-436 

-1391 

-1397 

-1385 

-1382 

-439 

-1389 

-1396 

-1393 

-1380 

5000 

-1724 

-550 

-1742 

-1748 

-1731 

-1727 

-551 

-1737 

-1747 

-1743 

-1723 

6000 

-2064 

-664 

-2093 

-2101 

-2074 

-2069 

-664 

-2083 

-2099 

-2091 

-2063 

7000 

-2399 

-781 

-2447 

-2458 

-2414 

•2404 

-780 

-2428 

-2455 

-2442 

-2399 

8000 

-2721 

-900 

-2805 

-2825 

-2749 

-2726 

-899 

-2771 

-2821 

-2797 

-2722 

8500 

•2868 

-963 

-2986 

-3017 

-2909 

-2868 

-964 

-2936 

-3011 

-2976 

-2869 

8600 

-2897 

-979 

-3027 

-3062 

-2943 

-2895 

-979 

-2971 

-3056 

-3016 

-2897 

8700 

-2922 

-993 

-3064 

-3105 

-2973 

-2917 

-994 

-3003 

-3098 

-3054 

-2922 

8800 

-2946 

-1008 

-3106 

-3151 

-3004 

-2940 

-1010 

-3038 

-3142 

-3094 

-2946 

8900 

-2973 

-1025 

-3151 

-3203 

-3038 

-2963 

-1029 

-3074 

-3193 

-3138 

-2972 

9000 

-2995 

-1042 

-3193 

-3253 

-3069 

-2981 

-1047 

-3106 

-3242 

-3180 

-2993 

NOTES:  1.  The  housing  assembly  consists  of  two  94  percent  alumina  ceramic  cylinders  with 

11.89  00  x  11.18  10  x  15  l  dimensions  bonded  to  a  simple  titaniun  joint  stiffener  Type  RJ1. 

2.  The  ends  of  the  ceramic  cylinders  were  closed  off  with  titanium  hemispheres  providing 
radial  and  axial  support. 

3.  The  assembly  successfully  withstood  11  pressure  cycles  to  9000  psi . 


60 


FEATURED  RESEARCH 


Table  2.  Strains  recorded  during  pressure  testing  of  case  3  test  assembly,  Sheet  2. 


|< . . . . Cylinder  No.  2  Midbay  Interior . . - . -->| 


Pressure 

Hoop 

Axial 

Hoop 

Hoop 

Hoop 

Hoop 

Axial 

Hoop 

Hoop 

Hoop 

Hoop 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1000 

-339 

-105 

-368 

-383 

-344 

-321 

-93 

-357 

-386 

-356 

-330 

2000 

-684 

-218 

-711 

-732 

-689 

-663 

-205 

-701 

-733 

-700 

-673 

3000 

-1032 

-331 

-1062 

-1087 

-1042 

-1009 

-320 

-1053 

-1088 

-1051 

-1021 

4000 

-1376 

-444 

-1405 

-1437 

-1387 

-1347 

-434 

-1396 

-1437 

-1394 

-1362 

5000 

-1723 

-558 

-1755 

-1793 

-1737 

-1690 

-549 

-1745 

-1793 

-1742 

-1707 

6000 

-2068 

-673 

-2103 

-2151 

-2085 

-2027 

-667 

-2090 

-2151 

-2088 

-2049 

7000 

-2413 

-793 

-2452 

-2514 

-2432 

-2371 

-788 

-2436 

-2516 

-2434 

-2386 

8000 

-2750 

-917 

-2804 

-2887 

-2776 

-2703 

-918 

-2779 

-2893 

-2779 

-2710 

8500 

-2901 

-976 

-2982 

-3081 

-2941 

-2846 

-983 

-2943 

-3091 

-2952 

-2855 

8600 

-2984 

•1045 

-3021 

-3127 

-2976 

-2925 

-1052 

-2977 

-3138 

-2991 

-2883 

8700 

-2957 

-1004 

-3058 

-3170 

-3008 

-2899 

-1014 

-3009 

-3182 

-3027 

-2907 

8800 

-2985 

-1021 

-3098 

-3216 

-3040 

-2901 

-1032 

-3041 

-3229 

-3064 

-2931 

8900 

-3014 

-1055 

-3142 

-3267 

-3076 

-2918 

-1054 

-3078 

-3283 

-3106 

-2957 

9000 

-3043 

•1089 

-3184 

-3317 

•3108 

-2935 

-1076 

-3109 

-3335 

-3146 

-2977 
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Table  2.  Strains  recorded  during  pressure  testing  of  case  3  test  assembly.  Sheet  3. 


l< . 

Cyl  i  rider 

No.  1  & 

2  Midbay 

Exterior  - 

. >1 

Pressure 

Hoop 

Axial 

Hoop 

Axial 

Hoop 

Axial 

Hoop 

Axial 

34 

35 

36 

37 

38 

39 

40 

41 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1000 

-322 

-93 

-319 

-99 

-337 

-102 

-361 

-83 

2000 

-634 

-198 

-630 

-204 

-646 

-205 

-675 

-193 

3000 

-953 

-305 

-947 

-308 

-966 

-314 

-999 

-305 

4000 

-1265 

-410 

-1258 

-413 

-1280 

-421 

-1317 

-414 

5000 

-1585 

-517 

-1573 

-517 

-1600 

-531 

-1642 

-527 

6000 

-1905 

-627 

-1890 

-625 

-1918 

-642 

-1968 

-642 

7000 

-2230 

-742 

-2213 

-738 

-2240 

-757 

-2301 

-760 

8000 

-2562 

•866 

-2548 

-861 

-2570 

-877 

-2649 

-889 

8500 

-2732 

-926 

-2718 

-923 

-2730 

-937 

-2826 

-960 

8600 

-2772 

-946 

-2813 

-939 

-2814 

-950 

-2866 

-979 

8700 

-2809 

-958 

-2839 

-956 

-2874 

-968 

-2909 

-996 

8800 

-2848 

-976 

-2841 

-974 

-2844 

-985 

-2954 

-1015 

8900 

-2893 

-998 

-2891 

-994 

-2887 

-1005 

-3006 

-1037 

9000 

-2945 

-1033 

-2941 

-1014 

-2930 

-1025 

-3058 

-1059 
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Table  2.  Strains  recorded  during  pressure  testing  of  case  3  test  assembly,  Sheet  4. 


l< . 

■--  Joint 

Stiffener, 

Internal 

Flange  Center  — 

. >1 

Pressure 

Hoop 

Axial 

Hoop 

Hoop 

Hoop 

Hoop 

Axial 

Hoop 

Hoop 

Hoop 

Hoop 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1000 

-176 

201 

-224 

-226 

-174 

-143 

204 

-197 

-255 

-215 

-170 

2000 

-369 

431 

-436 

-443 

-384 

-330 

447 

-409 

-486 

-426 

-372 

3000 

-570 

674 

-658 

-670 

-604 

-524 

698 

-625 

-720 

-644 

-579 

4000 

-765 

908 

-875 

-897 

-815 

-708 

935 

•830 

-946 

-859 

-775 

5000 

-962 

1146 

-1097 

-1131 

-1024 

-886 

1162 

-1027 

-1171 

-1081 

-964 

6000 

-1154 

1380 

-1323 

-1377 

-1239 

-891 

1383 

-1232 

-1416 

-1302 

-1139 

7000 

-1339 

1610 

-1568 

-1655 

-1460 

-1060 

1591 

-1441 

-1691 

-1531 

-1302 

8000 

-1493 

1807 

-1846 

-1986 

-1668 

-1149 

1752 

-1644 

-2026 

-1788 

-1421 

8500 

-1536 

1873 

-2014 

-2203 

-1757 

-1141 

1776 

-1733 

-2248 

-1938 

-1426 

8600 

-1532 

1876 

-2059 

-2267 

-1775 

-1124 

1763 

-1746 

-2312 

-1977 

-1410 

8700 

-1527 

1876 

-2104 

-2332 

-1791 

-1093 

1748 

-1760 

-2377 

-2017 

-1392 

8800 

-1521 

1876 

-2155 

-2405 

-1806 

-1057 

1728 

-1775 

-2452 

-2061 

-1368 

8900 

-1510 

1871 

-2217 

-2493 

-1820 

-1026 

1700 

-1791 

-2543 

-2114 

-1335 

9000 

-1495 

1862 

-2278 

-2584 

-1830 

-969 

1665 

-1805 

-2637 

-2167 

-1295 
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Table  2.  Strains  recorded  during  pressure  testing  of  case  3  test  assembly,  Sheet  5. 


| < -  Joint  Stiffener,  Internal  Flange  Edges - >| 


Pressure 

Hoop 

Axial 

Hoop 

Axial 

Hoop 

Axial 

Hoop 

Axial 

42 

43 

44 

45 

46 

47 

48 

49 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1000 

-211 

66 

-90 

32 

-142 

55 

-101 

36 

2000 

-415 

134 

-220 

80 

-312 

125 

-242 

90 

3000 

-610 

200 

-363 

134 

-486 

194 

-387 

147 

4000 

-795 

263 

-500 

184 

-649 

260 

-529 

202 

5000 

-975 

326 

-638 

236 

-668 

317 

-682 

262 

6000 

-1144 

385 

-776 

286 

-807 

371 

-822 

317 

7000 

-1291 

437 

-911 

337 

-902 

422 

-945 

365 

8000 

-1393 

476 

-1014 

374 

-980 

458 

-1017 

395 

8500 

-1401 

482 

-1029 

380 

-967 

456 

-1001 

392 

8600 

-1384 

478 

-1021 

378 

-941 

448 

-984 

386 

8700 

-1367 

473 

-1012 

374 

-916 

441 

-963 

379 

8800 

-1348 

468 

-999 

370 

-889 

431 

-936 

371 

8900 

-1323 

461 

-980 

363 

-850 

418 

-901 

359 

9000 

-1294 

453 

-956 

354 

-807 

404 

-858 

344 
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Table  2.  Strains  recorded  during  pressure  testing  of  case  3  test  assembly,  Sheet  6. 


Joint  Stiffener,  External  Flange 


\< . 

--  Center 

. >l<- 

--  Edges 

. *1 

Pressure 

Hoop 

Axial 

Hoop 

Axial 

Hoop 

Axial 

Hoop 

Axial 

Hoop 

Axial 

Hoop 

Axial 

SO 

51 

52 

S3 

54 

55 

56 

57 

58 

59 

60 

61 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1000 

-139 

38 

-188 

138 

-241 

-1 

-236 

0 

-238 

-9 

-224 

34 

2000 

-277 

73 

-347 

221 

-454 

-1 

-439 

0 

-454 

-10 

-419 

70 

3000 

-425 

130 

-512 

310 

-671 

0 

-645 

0 

-675 

-5 

-62 2 

109 

4000 

-571 

185 

-668 

414 

-887 

1 

-842 

0 

-890 

5 

-824 

143 

5000 

-727 

255 

-815 

539 

-1111 

3 

-1042 

0 

-1105 

34 

-1038 

170 

6000 

-888 

352 

-972 

659 

-1333 

6 

-1243 

0 

-1326 

66 

-1257 

205 

7000 

-1058 

485 

-1144 

780 

-1563 

10 

-1450 

0 

-1565 

107 

-1490 

249 

8000 

-1242 

602 

-1351 

903 

-1813 

16 

-1683 

0 

-1834 

154 

-1751 

307 

8500 

-1349 

669 

-1483 

973 

-1958 

21 

-1817 

0 

-1997 

190 

-1910 

342 

8600 

-1381 

704 

-1517 

1008 

•2000 

23 

-1851 

0 

-2044 

209 

-1957 

351 

8700 

-1410 

735 

-1555 

1039 

-2040 

2S 

-1884 

0 

-2089 

226 

-2001 

362 

8800 

-1443 

762 

-1599 

1070 

•2081 

27 

-1920 

0 

-2140 

243 

-2050 

375 

8900 

-1482 

788 

-1652 

1101 

-2128 

28 

-1965 

0 

-2198 

262 

•2108 

390 

9000 

-1519 

810 

-1708 

1128 

-2175 

30 

-2008 

0 

-2257 

280 

-2167 

405 
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APPENDIX  A:  NRaD  EPOXY  BONDING 
PROCEDURES  FOR  METALLIC  JOINT 
RINGS  USED  IN  CERAMIC  UNDER¬ 
WATER  PRESSURE  HOUSINGS 


A-1 
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FIGURE 

A-1 .  Axial  clearance  spacer  for  epoxy  bonding. 
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APPENDIX  A:  NRaD  EPOXY  BONDING 
PROCEDURES  FOR  METALLIC  JOINT 
RINGS  USED  IN  CERAMIC  UNDER¬ 
WATER  PRESSURE  HOUSINGS 


1 .  Wipe  bond  surfaces  of  the  ceramic  hull  com¬ 
ponent  with  a  clean  cloth  using  methyl  ethyl 
ketone  (MEK)  until  doth  shows  no  further  dis¬ 
coloration. 

2.  Wipe  the  interior  surfaces  of  the  metallic  joint 
ring  with  a  clean  cloth  using  MEK  until  cloth 
shows  no  further  discoloration. 

3.  For  titanium  joint  rings:  Passivate  the  interior 
surfaces  of  the  joint  ring  by  applying  a  layer  of 
PASA  JELL  107  and  allowing  it  to  etch  the 
titanium  surfaces  for  30  minutes.  Thoroughly 
rinse  off  interior  of  the  metallic  end  cap  and 
allow  surfaces  of  titanium  to  air  dry.  Air ’"g 
can  be  accelerated  with  a  forced-air  heater. 

4.  Mix  1 00  parts  CIBA  Geigy  601 0  epoxy  resin 
with  70  parts  CIBA  Geigy  283  hardener.  Use 
a  vacuum  chamber  to  remove  air  bubbles 
introduced  to  epoxy  during  mixing. 

5.  Pour  a .  1 2-inch-thick  layer  of  the  epoxy  mix¬ 
ture  into  the  interior  of  the  cylinder  end  cap. 
Place  axial  dearance  spacer  on  top  of  the 
epoxy  layer  and  press  it  down  through  the 
epoxy  mixture  to  the  bottom  of  the  joint  ring 
using  a  clean  tool.  A  typical  NRaD  spacer  giv¬ 
ing  an  axial  clearance  of  .010  of  an  inch  for  a 
12-inch  OD  ceramic  cylinder  (.412-inch-thick 
wall)  or  hemisphere  is  shown  in  figure  A-1 . 


6.  Pour  additional  epoxy  mixture  into  the  bottom 
of  the  joint  ring  until  the  epoxy  fills  at  least  half 
the  depth  of  the  interior  of  the  joint  ring. 

7.  Lower  the  end  of  the  ceramic  hull  component 
into  the  joint  ring  interior  (partially  filled  with 
the  epoxy  mixture)  keeping  the  ceramic  hull 
centered  within  the  joint  ring.  Allow  the  hull  to 
settle  evenly  into  the  joint  ring  interior  until  the 
hull  comes  to  rest  on  the  spacer  at  the  bottom 
of  the  joint  ring  interior.  Additional  weight  can 
be  placed  on  top  of  ceramic  hull  component  to 
help  it  settle  evenly  through  the  epoxy  mix¬ 
ture.  Care  should  be  taken  to  assure  the  hull 
remans  centered  with  the  joint  ring  and  the 
center  line  of  the  hull  remains  perpendicular  to 
the  working  surface  throughout  the  bonding 
procedure. 

8.  Leave  the  bonded  assembly  and  any  settling 
weight  used  undisturbed  for  at  least  24  hours 
to  allow  the  epoxy  mixture  to  begin  to  cure. 
Once  this  initial  24-hour  period  has  passed, 
completely  remove  any  excess  epoxy  mixture 
remaining  on  the  exterior  surface  of  the  joint 
ring  that  extruded  out  during  assembly. 
Application  of  a  thin  coat  of  silicone-based 
mold  release  to  the  exterior  surfaces  of  the 
joint  ring  prior  to  bonding  the  hull  may  be  used 
to  aid  in  the  subsequent  cleanup.  Care  should 
be  taken  to  not  allow  mold  release  on  or  near 
any  bonding  surfaces  on  the  hull  or  joint  ring. 

9.  After  the  epoxy  has  cured,  apply  a  bead  of 
R7V  or  other  suitable  sealant  to  externally 
exposed  portions  of  the  epoxy  bond  between 
the  metallic  joint  ring  and  outer  surface  of  the 
ceramic  hull  component. 


A-3 


Figure  A-1 .  Axial  clearance  spacer  for  epoxy  bonding. 
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APPENDIX  B:  B0S0R4  BUCKLING 
ANALYSIS  INPUT  FILE  FOR  CASE  4 
TEST  ASSEMBLY 


B-1 


FEATURED  RESEARCH 


FIGURES 

B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
3-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 
B-1 .  Scale  model  housing 


data  for  case  4,  Sheet  1 . 
data  for  case  4,  Sheet  2. 
data  for  case  4,  Sheet  3. 
data  for  case  4,  Sheet  4. 
data  for  case  4,  Sheet  5. 
data  for  case  4,  Sheet  6. 
data  for  case  4,  Sheet  7. 
data  for  case  4,  Sheet  8. 
data  for  case  4,  Sheet  9. 
data  for  case  4,  Sheet  10. 
data  for  case  4,  Sheet  1 1 . 
data  for  case  4,  Sheet  12. 
data  for  case  4,  Sheet  13. 
data  for  case  4,  Sheet  14. 
data  for  case  4,  Sheet  15. 


B-2 


FEATURED  RESEARCH 


CASE  4 :  SCALE  MODEL  HOUSING 

1  $  INDIC  -  analysis  type  indicator 

1  $  NPRT  -  output  options  ( 1 -minimum,  2-medium,  3-maximum) 

0  S  ISTRES-  output  control  (O-resultants,  1-sigma,  2-epsilon) 

15  $  NSEG  -  number  of  shell  segments  (less  than  95) 

H  $ 

H  S  SEGMENT  NUMBER  11111111 

H  $  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

10  $  NMESH  -  number  of  node  points  (5  -  min.;  98  -  max.)(  1) 

3  $  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 

H  $  REFERENCE  SURFACE  GEOMETRY  FOLLOWS. . . 

2  S  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 

0.  5  Rl  -  radius  at  beginning  of  segment  (see  p.  66) 

0.  $  zi  -  global  axial  coordinate  at  beginning  of  segment 

5.353000  S  R2  -  radius  at  end  of  segment 

3.30 6000  $  Z2  -  global  axial  coordinate  at  end  of  segment 

0.  $  RC  -  radius  from  axis  of  rev.  to  center  of  curvature 

5 . 988000  $  ZC  -  axial  coordinate  of  center  of  curvature 

-1  $  SROT-indicator  for  direction  of  increasing  arc  (-1.  or  +1.) 

H  $  IMPERFECTION  SHAPE  FOLLOWS... 

0  $  IMP  -  indicator  for  imperfection  (0-none,  1-some) 

H  $  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 

3  $  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 

0.2710000  $  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

N  $  Do  you  want  to  print  out  r(s),  r’(s),  etc.  for  this  segment? 

H  $  DISCRETE  RING  INPUT  FOLLOWS... 

0  5  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 

0  $  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

H  $  LINE  LOAD  INPUT  FOLLOWS . . . 

0  $  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 

H  $  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

1  $  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 

H  $  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 

1  $  NLTYPE- control  (0,1, 2, 3)  for  type  of  surface  loading 

2  $  NPSTAT-  number  of  meridional  callouts  for  surface  loading 

0  $  NLOAD(l) -indicator  for  meridional  traction  (0-none,  l-some) 

0  S  NI/3AD( 2) -indicator  for  circumferential  traction 

1  $  NLOAD(3 )-indicator  for  normal  pressure  (0-none,  1-some) 

-1  $  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

-1  $  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  2) 

2  5  NTYPE  -  control  for  meaning  of  loading  callout  (2-z,  3»r) 

0.  5  Z( I )  -  axial  coordinate  of  Ith  loading  callout,  z(  1) 

3.306000  5  Z(I)  -  axial  coordinate  of  Ith  loading  callout,  z(  2) 

H  S  SHELL  WALL  CONSTRUCTION  FOLLOWS... 

2  5  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 

0. 4400000E+08  $  E  -  Young's  modulus  for  skin 

0.2100000  $  v  -  Poisson's  ratio  for  skin 

0.  S  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec**2/in**4 ) 

0.  $  ALPHA  -  coefficient  of  thermal  expansion 

0  $  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 

1  $  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

N  $  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

N  $  Do  you  want  to  print  out  distributed  loads  along  meridian? 

H  $ 

H  $  SEGMENT  NUMBER  2  2222222 

H  $  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

5  $  NMESH  -  number  of  node  points  (5  -  min.;  98  -  max.)(  2) 

3  $  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 

H  S  REFERENCE  SURFACE  GEOMETRY  FOLLOWS... 

1  S  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 

5.353000  $  Rl  -  radius  at  beginning  of  segment  (see  p.  66) 

3.306000  $  Zl  -  global  axial  coordinate  at  beginning  of  segment 

6.024000  S  R2  -  radius  at  end  of  segment 

4.783000  S  12  -  global  axial  coordinate  at  end  of  segment 

H  $  IMPERFECTION  SHAPE  FOLLOWS... 

C  S  IMP  »  indicator  for  imperfection  (0-none,  1-some) 


B-1 .  Scale  model  housing  data  for  case  4,  Sheet  1 . 
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S  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 
$  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 
$  NZVALU-  number  of  meridional  callouts  for  ref.  surf. 

$  NTYPE  -  control  for  meaning  of  callout  ( 2-z ,  3-r) 

$  Z(I)  -  axial  coordinate  of  Ith  callout,  z(  1) 

$  Z(I)  =  axial  coordinate  of  Ith  callout,  z(  2) 

$  ZVAL  -  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  i) 

$  ZVAL  -  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  2) 

$  Do  you  want  to  print  out  r(s),  r’(s),  etc.  for  this  segment? 
$  DISCRETE  RING  INPUT  FOLLOWS... 

$  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 
$  K-elastic  foundation  modulus  (e.g.  lb/in*»3)in  this  seg. 

S  LINE  LOAD  INPUT  FOLLOWS... 

$  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 
$  DISTRIBUTED  LOAD  INPUT  FOLLOWS . . . 

$  IDISAB-  indicator  ( 0 ,  1 ,  2  or  3 )  for  load  set  A  and  B 
S  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 
$  NLTYPE-control  (0,1, 2, 3)  for  type  of  surface  loading 
$  NPSTAT-  number  of  meridional  callouts  for  surface  loading 
$  NLOAD(l ^indicator  for  meridional  traction  (0-none,  1-some) 

$  NL0AD( 2) -indicator  for  circumferential  traction 
S  NLOAD( 3) -indicator  for  normal  pressure  (0-none,  1-some) 

5  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

$  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  2) 

S  NTYPE  -  control  for  meaning  of  loading  callout  (2-z,  3-r) 

$  Z(I)  -  axial  coordinate  of  Ith  loading  callout,  z(  1) 

$  Z(I)  -  axial  coordinate  of  Ith  loading  callout,  z(  2) 

S  SHELL  WALL  CONSTRUCTION  FOLLOWS... 

$  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 
S  E  -  Young ' s  modulus  for  shin 

$  U  -  Poisson's  ratio  for  shin 

$  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec ** 2/in » *4 ) 
$  ALPHA  -  coefficient  of  thermal  expansion 
S  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 
S  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

$  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

5  Do  you  want  to  print  out  distributed  loads  along  meridian? 

S 

S  SEGMENT  NUMBER  33333333 
S  NODAL  POINT  DISTRIBUTION  FOLLOWS. . . 

$  NMESH  -  number  of  node  points  (5  -  min.;  98  -  max.)(  3) 

S  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  SDacing 
$  REFERENCE  SURFACE  GEOMETRY  FOLLOWS. 

S  NS H ape-  indicator  (1,2  or  4)  for  geometry  of  meridian 
S  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

S  Z1  -  global  axial  coordinate  at  beginning  of  segment 

S  R2  -  radius  at  end  of  segment 

$  Z2  «  global  axial  coordinate  at  end  of  segment 

S  IMPERFECTION  SHAPE  FOLLOWS... 

S  IMP  -  indicator  for  imperfection  (0-none,  1-some) 

$  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 
S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 
S  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

S  Do  you  want  to  print  out  r(s),  r’(s),  etc.  for  this  segment? 
S  DISCRETE  RING  INPUT  FOLLOWS... 

S  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 
S  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

S  LINE  LOAD  INPUT  FOLLOWS... 

$  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 
$  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

S  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
$  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 
$  NLTYPE-control  (0,1, 2, 3)  for  type  of  surface  loading 
$  NPSTAT-  number  of  meridional  callouts  for  surface  loading 
$  NLCAD;  1) -indicator  for  meridional  traction  (0-none,  1-soir.e) 

$  NL0AD( 2 ) -indicator  for  circumferential  traction 


B-1 .  Scale  model  housing  data  for  case  4,  Sheet  2. 
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$  NLOAD(3)-indicator  for  normal  pressure  (0-none,  1-some) 

$  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

S  PN(i)  -  normal  pressure  ( p.74)  at  ith  callout,  PN(  2) 

5  NTYPE  -  control  for  meaning  of  loading  callout  (2-z,  3«r) 

$  Z( I )  =  axial  coordinate  of  Ith  loading  callout,  z(  1) 

$  Z ( I )  -  axial  coordinate  of  Ith  loading  callout,  z(  2) 

S  SHELL  WALL  CONSTRUCTION  FOLLOWS ... 

?  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 
$  E  -  Young's  modulus  for  skin 

S  U  -  Poisson '  s  ratio  for  skin 

$  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec**2/in**4 ) 
S  ALPHA  -  coefficient  of  thermal  expansion 
$  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 

$  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

$  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

$  Do  you  want  to  print  out  distributed  loads  along  meridian? 

S 

S  SEGMENT  NUMBER  44444444 
$  NODAL  POINT  DISTRIBUTION  FOLLOWS. . . 

$  NMESH  =  number  of  node  points  (5  -  min.,-  98  -  rnax.M  4) 

$  NTYPEH-  control  integer  (1  or  3 )  for  nodal  point  spacing 
$  REFERENCE  SURFACE  GEOMETRY  FOLLOWS ... 

$  NSHAPE-  indicator  <1,2  or  4)  for  geometry  of  meridian 
$  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

$  Z1  -  global  axial  coordinate  at  beginning  of  segment 

$  R2  -  radius  at  end  of  segment 

$  Z2  -  global  axial  coordinate  at  end  of  segment 

S  IMPERFECTION  SHAPE  FOLLOWS ... 

5  IMP  -  indicator  for  imperfection  (0=none,  1-some) 

S  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 
$  NTYFEZ-  control  (1  or  3)  for  reference  surface  location 
S  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

S  Do  you  want  to  print  out  r{s),  r'(s),  etc.  for  this  segment? 
$  DISCRETE  RING  INPUT  FOLLOWS ... 

$  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 
$  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

$  LINE  LOAD  INPUT  FOLLOWS. . . 

S  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 
S  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

S  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
$  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 
S  NLTYPE-control  (0,1, 2, 3)  for  type  of  surface  loading 
S  NPSTAT-  number  of  meridional  callouts  for  surface  loading 
S  NLOAD(l) -indicator  for  meridional  traction  (0-none,  1-some) 

5  NL0AD(2)-indicator  for  circumferential  traction 
S  NL0AD( 3 )-indicator  for  normal  pressure  (0-none,  1-some) 

$  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

S  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  2) 

S  NTYPE  -  control  for  meaning  of  loading  callout  (2-z,  3-r) 

S  Z(I)  -  axial  coordinate  of  Ith  loading  callout,  z(  1) 

$  Z(I)  -  axial  coordinate  of  Ith  loading  callout,  z(  2) 

S  SHELL  WALL  CONSTRUCTION  FOLLOWS... 

$  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 
S  E  -  Young ' s  modulus  for  skin 

S  U  -  Poisson's  ratio  for  skin 

S  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb- sec** 2/in ’-*4  ) 
S  ALPHA  -  coefficient  of  thermal  expansion 
S  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 
$  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

S  Do  you  want  to  print  out  ref.  surf,  location  and  thickness’ 

S  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations’ 

$  Do  you  want  to  print  out  distributed  loads  along  meridian? 

$ 

$ 

$  SEGMENT  NUMBER  55555555 
$  NODAL  POINT  DISTRIBUTION  FOLLOWS. . . 


B-1 .  Scale  model  housing  data  for  case  4,  Sheet  3. 
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13  $  NMESH  -  number  of  node  points  (5  -  min.;  98  -  max.)(  5) 

3  S  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 

H  $  REFERENCE  SURFACE  GEOMETRY  FOLLOWS... 

1  S  NSHAPE=  indicator  (1,2  or  4)  for  geometry  of  meridian 

5.804000  $  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

21.57300  $  Z1  -  global  axial  coordinate  at  beginning  of  segment 

5.804000  S  R2  -  radius  at  end  of  segment 

21 . 81300  $  Z2  *  global  axial  coordinate  at  end  of  segment 

H  $  IMPERFECTION  SHAPE  FOLLOWS... 

0  $  IMP  -  indicator  for  imperfection  (0=none,  1-some) 

H  $  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 

3  S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 

0.2340000  $  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

N  $  Do  you  want  to  print  out  r(s),  r'(s),  etc.  for  this  segment’ 

H  $  DISCRETE  RING  INPUT  FOLLOWS ... 

0  S  NRINGS-  number  (max* 20)  of  discrete  rings  in  this  segment 

0  $  K=elastic  foundation  modulus  (e.g.  lb/in**3 . in  this  seg. 

H  S  LINE  LOAD  INPUT  FOLLOWS ... 

0  $  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 

H  $  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

1  $  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 

H  $  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 

1  $  NLTYPE-control  (0,1, 2, 3)  for  type  of  surface  loading 

2  $  NPSTAT-  number  of  meridional  callouts  for  surface  loading 

0  S  NLOAD(l)=indicator  for  meridional  traction  (0-none,  l-sctne) 

0  $  NLOAD(2)=indicator  for  circumferential  traction 

1  $  NLOAD(3)-indicator  for  normal  pressure  (0=none,  l=some) 

-1  S  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

-1  $  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  2) 

2  S  NTYPE  -  control  for  meaning  of  loading  callout  (2-z,  3-r) 

21.57300  $  Z(I)  =  axial  coordinate  of  Ith  loading  callout,  z(  1) 

21.81300  $  Z( I )  -  axial  coordinate  of  Ith  loading  callout,  z(  2) 

H  $  SHELL  WALL  CONSTRUCTION  FOLLOWS ... 

2  $  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 

0 . 1640000E+08  5  E  -  Young's  modulus  for  skin 

0.3100000  $  U  -  Poisson's  ratio  for  skin 

0.  $  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec** 2/in* *4 ) 

0 .  S  ALPHA  -  coefficient  of  thermal  expansion 

0  $  NRS  =  control  (0  or  1)  for  addition  of  smeared  stiffeners 

0  S  NSUR  -  control  for  thickness  input  (0  or  1  or  -1 ) 

Y  5  Do  you  want  to  print  out  ref.  surf,  location  and  thickness? 

N  5  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

N  S  Do  you  want  to  print  out  distributed  loads  along  meridian? 

H  S 

H  S  SEGMENT  NUMBER  66666666 

H  S  NODAL  POINT  DISTRIBUTION  FOLLOWS. . . 

13  $  NMESH  -  number  of  node  points  (5  -  min.;  98  -  max. ) (  6) 

3  S  NTYPEH—  control  integer  (1  or  3)  for  nodal  point  spacing 

H  $  REFERENCE  SURFACE  GEOMETRY  FOLLOWS ... 

1  5  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 

5.510000  S  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

20^61000  S  Z1  =  global  axial  coordinate  at  beginning  of  segment 

5.510000  S  R2  -  radius  at  end  of  segment 

22.77600  S  Z2  -  global  axial  coordinate  at  end  of  segment 

H  $  IMPERFECTION  SHAPE  FOLLOWS... 

0  5  IMP  -  indicator  for  imperfection  (0-none,  1-some; 

H  $  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 

1  S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 

4  $  NZVALU-  number  of  meridional  callouts  for  ref.  surf. 

2  $  NTYPE  -  control  for  meaning  of  callout  (2-z,  3-r) 

20.61000  ?  Z(I)  -  axial  coordinate  of  Ith  callout,  z(  1) 

21.57300  $  Z(I)  -  axial  coordinate  of  Ith  callout,  z(  2) 

21.81300  S  2(1)  -  axial  coordinate  of  Ith  callout,  z(  3) 

22.77600  S  2(1)  -  axial  coordinate  of  Ith  callout,  z(  4) 

0.  S  ZVAL  -  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  1) 

0 . 6000000E-01  S  ZVAL  =  distance  from,  leftmost  surf,  to  ref.  surf.,ZVAL(  2) 

B-1.  Scale  model  housing  data  for  case  4,  Sheet  4. 
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ZVAL  -  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  3) 
ZVAL  -  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  4) 
Do  you  want  to  print  out  r(s),  r'(s),  etc.  for  this  segment? 
DISCRETE  RING  INPUT  FOLLOWS... 

NRINGS-  number  (max- 20)  of  discrete  rings  in  this  segment 
K=elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

LINE  IOAD  INPUT  FOLLOWS.. . 

LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 
DISTRIBUTED  LOAD  INPUT  FOLLOWS. . . 

IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
SHELL  WALL  CONSTRUCTION  FOLLOWS. . . 

NWALL- index  (1,  2,  4,  5,  6,  1,  8)  for  wall  construction 
E  -  Young's  modulus  for  skin 

U  -  Poisson  1 s  ratio  for  skin 

SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec-*2/in»-4 ) 
ALPHA  -  coefficient  of  thermal  expansion 

NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 
NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

NTYPET-  index  (1  or  3)  for  type  of  input  for  thickness 
NTVALU-  number  of  callouts  along  segment  for  thickness 
NTYPE  -  control  for  meaning  of  thickness  callout  (2-z,  3=r) 

Z ( I )  -  axial  coordinate  of  Ith  thickness  callout,  z(  1) 

Z(I)  -  axial  coordinate  of  Ith  thickness  callout,  z(  2) 

Z(I)  -  axial  coordinate  of  Ith  thickness  callout,  z(  3) 

Z < I )  -  axial  coordinate  of  Ith  thickness  callout,  z(  4) 

TVAL(i)  -  thickness  at  Ith  callout,  TVAL(  1) 

TVAL(i)  -  thickness  at  Ith  callout,  TVALC  2) 

TVAL(i)  -  thickness  at  Ith  callout,  TVAL(  3) 

TVAL(i)  -  thickness  at  Ith  callout,  TVAL(  4) 

Do  you  want  to  print  out  ref.  surf,  location  and  thickness'’ 
Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations’ 
Do  you  want  to  print  out  distributed  loads  along  meridian’ 


H  S  SEGMENT  NUMBER  7’77’777 

r  S  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

9  S  NXESH  -  number  of  node  points  (5  -  min.,-  98  -  max.)(  7) 

i  s  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 

K  S  REFERENCE  SURFACE  GEOMETRY  FOLLOWS . . . 

1  S  '.'SHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 

5.4  531’C1  S  Rl  -  radius  at  beginning  of  segment  (see  p.  66; 

21.693CC  S  21  -  global  axial  coordinate  at  beginning  of  segment 

303302  S  R2  -  radius  at  end  of  segment 

21 .69300  S  22  -  global  axial  coordinate  at  end  of  segment 

H  S  IMPERFECTION  SHAPE  FOLLOWS... 

C  S  IX?  »  indicator  for  imperfection  (C-none,  i-some) 

K  5  REFERENCE  SURFACE  LOCATION  RELATIVE  10  WALL 

:  S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 

C  .  12CL0  2C  S  ZVAL  -  distance  from,  leftmost  surf,  to  reference  surf. 

S  Do  you  want  to  print  out  r(s),  r'(s),  etc.  for  this  seament’ 
H  $  DISCRETE  RING  INPUT  FOLLOWS. . . 

'  S  NRINGS-  number  (max-2C)  of  discrete  rings  m  this  segment 

C  S  K-eiastic  foundation  modulus  (e.c.  lb/in**3)in  this  seg. 

K  S  LINE  LOAD  INPUT  FOLLOWS... 

2  $  LINTYP-  indicator  (0,  1,  2  or  3;  for  tvpe  of  line  loads 

K  5  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

C>  S  IDISAB-  indicator  (0.  1,  2  or  3)  fcr  load  set  A  and  B 

K  $  SHELL  WALL  CONSTRUCTION  FOLLOWS... 

2  $  NWALL-mdex  (1.  2,  4.  5,  6,  ",  3)  for  wall  construction 

3. 164C03CE-08  S  E  -  Young's  modulus  for  skin 

3.3100332  $  U  -  Poisson's  ratio  for  skin 

3.  S  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec* *2/m**4 ) 

3  S  ALPHA  -  coefficient  of  thermal  expansion 

2  S  NRS  -  control  ;3  or  1)  for  addition  of  smeared  stiffeners 

2  S  NSUR  -  control  for  thickness  input  (C  or  1  or  -1: 

Y  5  Do  yci  want  to  print  out  ref  surf,  location  and  thickness’ 

S  ?  Do  you  want  to  print  out  the  C(i.;)  at  meridional  stations’ 


B-1 .  Scale  model  housing  data  for  case  4,  Sheet  5. 
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N  $  Do  you  want  to  print  out  distributed  loads  along  meridian? 

H  $ 

H  $  SEGMENT  NUMBER  88888888 

H  $  NODAL  POINT  DISTRIBUTION  FOLLOWS. . . 

13  S  NMESH  -  number  of  node  points  (5  =  min.;  98  -  max. )(  8) 

3  S  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 

H  $  REFERENCE  SURFACE  GEOMETRY  FOLLOWS... 

1  $  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 

5.105000  $  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

20.61000  $  Z1  -  global  axial  coordinate  at  beginning  of  segment 

5.105000  $  R2  -  radius  at  end  of  segment 

22.77600  S  Z2  =  global  axial  coordinate  at  end  of  segment 

H  S  IMPERFECTION  SHAPE  FOLLOWS... 

0  S  IMP  -  indicator  for  imperfection  (0-none,  1-some) 

H  $  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 

3  S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 

0.1650000  S  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

N  $  Do  you  want  to  print  out  r(s),  r'ts),  etc.  for  this  segment? 

H  S  DISCRETE  RING  INPUT  FOLLOWS... 

0  S  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 

0  $  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

H  5  LINE  LOAD  INPUT  FOLLOWS . . . 

0  S  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 

H  S  DISTRIBUTED  LOAD  INPUT  FOLLOWS . . . 

0  S  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 

H  S  SHELL  WALL  CONSTRUCTION  FOLLOWS. . . 

2  S  NWALL- index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 

0 . 1640000E^08  S  E  «  Young's  modulus  for  skin 

0.3100000  S  U  =  Poisson's  ratio  for  skin 

0.  S  SM  -mass  density  of  skin  (e.g.  alum. -. 00025  lb-sec»*2/in**4 ) 

0.  $  ALPHA  -  coefficient  of  thermal  expansion 

0  S  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 

-1  $  NSUR  -  control  for  thickness  input  (0  or  I  or  -1) 

1  S  NTYPET-  index  (1  or  3)  for  type  of  input  for  thickness 

4  S  NTVALU"  number  of  callouts  along  segment  for  thickness 

2  S  NTYPE  »  control  for  meaning  of  thickness  callout  (2«z,  3-r) 

20.61000  S  Z(I)  -  axial  coordinate  of  Ith  thickness  callout,  z(  1) 

21.57300  S  Z(I)  -  axial  coordinate  of  Ith  thickness  callout,  z(  2) 

21.81300  $  Z(I)  -  axial  coordinate  of  Ith  thickness  callout,  z(  3) 

22.77600  $  Z(I)  -  axial  coordinate  of  Ith  thickness  callout,  z(  4) 

0.3300000  5  TVAL(i)  -  thickness  at  Ith  callout,  TVAL(  1) 

0.3600000  $  TVAL(i)  -  thickness  at  Ith  callout,  TVAL(  2) 

0.3600000  5  TVAL(i)  -  thickness  at  Ith  callout,  TVAL(  3) 

0.3300000  S  TVAL(i)  -  thickness  at  Ith  callout,  TVAL(  4) 

Y  S  Dc  you  want  to  print  out  ref.  surf,  location  and  thickness? 

N  S  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

N  S  Do  you  want  to  print  out  distributed  loads  along  meridian? 

H  S 

H  S  SEGMENT  NUMBER  99999999 

H  S  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

5  S  NMESH  -  number  of  node  points  (5  -  min.,-  98  -  max.)(  9; 

3  S  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 

H  S  REFERENCE  SURFACE  GEOMETRY  FOLLOWS  ... 

1  S  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 

e.oeeooc  S  Si  -  radius  at  beginning  of  segment  (see  p.  66) 

20.61000  S  Z1  -  global  axial  coordinate  at  beginning  of  segment 

6.088000  S  R2  -  radius  at  end  of  segment 

20.85000  S  Z2  -  global  axial  coordinate  at  end  of  segment 

H  S  IMPERFECTION  SHAPE  FOLLOWS ... 

0  S  IMP  -  indicator  for  imperfection  (0-none,  1-some) 

H  S  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 

3  S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 

0 . 5000000E-0!  S  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

N  S  Do  you  want  to  print  out  r[s),  r'(s),  etc.  for  this  segment? 

H  S  DISCRETE  RING  INPUT  FOLLOWS... 

0  S  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 


B-1 .  ScaJe  model  housing  data  for  case  4,  Sheet  6. 
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$  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

$  LINE  LOAD  INPUT  FOLLOWS... 

?  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 
$  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

S  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
S  SHELL  WALL  CONSTRUCTION  FOLLOWS . . . 

$  NWALL- index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 
$  E  -  Young's  modulus  for  skin 

S  U  =  Poisson's  ratio  for  skin 

$  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec**2/in**4 ) 
S  ALPHA  -  coefficient  of  thermal  expansion 
$  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 
S  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

S  Do  you  want  to  print  out  ref.  surf,  location  and  thickness? 

S  Do  you  want  to  print  out  the  C(i, j )  at  meridional  stations? 

$  Do  you  want  to  print  out  distributed  loads  along  meridian’ 

$ 

S  SEGMENT  NUMBER  10  10  10  10  10  10  10  10 

S  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

$  NMESH  -  number  of  node  points  (5  -  min.;  98  -  max.) (10) 

S  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 
S  REFERENCE  SURFACE  GEOMETRY  FOLLOWS... 

$  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 
$  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

$  Z1  -  global  axial  coordinate  at  beginning  of  segment 

S  R2  -  radius  at  end  of  segment 

5  Z2  -  global  axial  coordinate  at  end  of  segment 

S  IMPERFECTION  SHAPE  FOLLOWS... 

S  IMP  -  indicator  for  imperfection  (0-none,  1-scme) 

S  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 
S  NTYPEZ*  control  (1  or  3 )  for  reference  surface  location 
S  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

S  Do  you  want  to  print  out  r(s),  r'(s),  etc.  for  this  segment? 
S  DISCRETE  RING  INPUT  FOLLOWS... 

S  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 
S  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

S  LINE  LOAD  INPUT  FOLLOWS. . . 

S  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 
S  DISTRIBUTED  LOAD  INPUT  FOLLOWS. . . 

S  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
S  SHELL  WALL  CONSTRUCTION  FOLLOWS... 

$  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 
S  E  -  Young ' s  modulus  for  skin 

$  U  -  Poisson's  ratio  for  skin 

S  SM  -mass  density  of  skin  (e.g.  alum.-. 00025  lb-sec** 2/in* *4) 
S  ALPHA  -  coefficient  of  thermal  expansion 
S  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 
S  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

S  Do  you  want  to  print  out  ref .  surf .  location  and  thickness? 

$  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

?  Do  you  want  to  print  out  distributed  loads  along  meridian? 

5 

S  SEGMENT  NUMBER  11  11  II  11  11  11  11  11 

S  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

S  NMESH  -  number  of  node  points  (5  -  min.,-  98  «  max.)  (11) 

S  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 
S  REFERENCE  SURFACE  GEOMETRY  FOLLOWS... 

S  NSHAPE-  indicator  (1,2  or  4 )  for  geometry  of  meridian 
S  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

S  Z1  -  global  axial  coordinate  at  beginning  of  segment 

S  R2  -  radius  at  end  of  segment 

S  Z2  -  global  axial  coordinate  at  end  of  segment 

S  IMPERFECTION  SHAPE  FOLLOWS... 

S  IMP  -  indicator  for  imperfection  (0-none,  1-some) 

S  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 
S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 


B-1 .  Scale  model  housing  data  for  case  4,  Sheet  7. 
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S  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

S  Do  you  want  to  print  out  r(s),  r’(s),  etc.  for  this  segment? 
$  DISCRETE  RING  INPUT  FOLLOWS... 

S  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 
$  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

S  LINE  LOAD  INPUT  FOLLOWS... 

S  LINTYP*  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 
S  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

S  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
S  SHELL  WALL  CONSTRUCTION  FOLLOWS . . . 

S  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 
$  E  =  young  ’  s '  modulus  for  skin 

$  U  -  Poisson's  ratio  for  skin 

$  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec**2/in**4 ) 
$  ALPHA  -  coefficient  of  thermal  expansion 

$  NRS  -  control  (0  or  l)  for  addition  of  smeared  stiffeners 
$  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

S  Do  you  want  to  print  out  ref.  surf,  location  and  thickness15 
S  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations5 
S  Do  you  want  to  print  out  distributed  loads  along  meridian? 

S  SEGMENT  NUMBER  12  12  12  12  12  12  12  12 

S  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

S  NMESH  -  number  of  node  points  (5  -  min.;  98  =  max.) (12) 

$  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 
S  REFERENCE  SURFACE  GEOMETRY  FOLLOWS ... 

$  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 
$  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

S  Z1  -  global  axial  coordinate  at  beginning  of  segment 

S  R2  -  radius  at  end  of  segment 

S  22  »  global  axial  coordinate  at  end  of  segment 

S  IMPERFECTION  SHAPE  FOLLOWS . . . 

$  IMP  =  indicator  for  imperfection  (0-none,  1-some) 

$  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 
S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 
S  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

$  Do  you  want  to  print  out  r(s),  r‘(s),  etc.  for  this  segment5 
$  DISCRETE  RING  INPUT  FOLLOWS... 

S  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 
S  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

S  LINE  LOAD  INPUT  FOLLOWS... 

S  LINTYP-  indicator  (0.  1,  2  or  3)  for  type  of  line  loads 
S  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

S  ID I SAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
S  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 
S  NLTYPE-control  ( 0 , 1 , 2 , 3 )  for  type  of  surface  loading 
S  NPSTAT-  number  of  meridional  callouts  for  surface  loading 
S  NLOAD(l) -indicator  for  meridional  traction  (0-none,  1-some) 

S  NLOAD; 2) -indicator  for  circumferential  traction 
S  NLOAD(3)-indicator  for  normal  pressure  (0-none,  1-some) 

S  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

S  PN'(i)  =  normal  pressure  (p.74)  at  ith  callout,  PN(  2) 

S  NTYPE  -  control  for  meaning  of  loading  callout  (2-z,  3-r) 

S  Z ( I )  -  axial  coordinate  of  Ith  loading  callout,  z(  1) 

S  Z(I)  -  axial  coordinate  of  Ith  loadina  callout,  z(  2) 

S  SHELL  WALL  CONSTRUCTION  FOLLOWS ... 

S  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 
S  E  -  Young's  modulus  for  skin 

S  U  -  Poisson's  ratio  for  skin 

S  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec**2/in**4 ) 
S  ALPHA  -  coefficient  of  thermal  expansion 
S  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 
5  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

S  Do  you  want  to  print  out  ref.  surf,  location  and  thickness5 
S  Dc  you  want  to  print  out  the  C(i,})  at  meridional  stations5 
S  Do  you  want  to  print  out  distributed  loads  along  meridian5 
$ 

B-1 .  Scale  model  housing  data  for  case  4,  Sheet  8. 
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H  $  SEGMENT  NUMBER  13  13  13  13  13  13  13  13 

H  $  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

5  $  NMESH  -  number  of  node  points  (S  «  min.;  98  -  max.) (13) 

3  5  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 

H  $  REFERENCE  SURFACE  GEOMETRY  FOLLOWS . . . 

1  S  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 

6.024000  s  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

37.39800  S  21  -  global  axial  coordinate  at  beginning  of  segment 

6.024000  S  R2  -  radius  at  end  of  segment 

38.60300  S  22  -  global  axial  coordinate  at  end  of  segment 

H  $  IMPERFECTION  SHAPE  FOLLOWS... 

0  $  IMP  -  indicator  for  imperfection  (0=none,  1-some) 

H  S  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 

3  $  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 

0.4340000  $  ZVAL  -  distance  from  leftmost  surf,  to  reference  surf. 

N  $  Do  you  want  to  print  out  r(s),  r'(s),  etc.  for  this  segment? 

H  $  DISCRETE  RING  INPUT  FOLLOWS . . . 

0  $  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 

0  $  K-elastic  foundation  modulus  (e.g.  lb/in*»3)in  this  seg. 

H  S  LINE  LOAD  INPUT  FOLLOWS ... 

0  5  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 

H  S  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

1  S  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 

H  $  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 

1  S  NLTYPE-control  (0,1, 2, 3)  for  type  of  surface  loading 

2  S  NPSTAT—  number  of  meridional  callouts  for  surface  loading 

0  S  NLOAD(l) -indicator  for  meridional  traction  (0-none,  1-some) 

0  S  nloadc 2) -indicator  for  circumferential  traction 

1  S  NL0AD( 3) -indicator  for  normal  pressure  (0-none,  1-some) 

-1  $  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

-1  $  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  2) 

2  $  NTYPE  -  control  for  meaning  of  loading  callout  (2-2,  3-r) 

37.39800  $  2(1)  -  axial  coordinate  of  Ith  loading  callout,  z(  1) 

38.60300  S  Z  ( I )  -  axial  coordinate  of  Ith  loading  callout,  z(  2) 

H  $  SHELL  WALL  CONSTRUCTION  FOLLOWS . . . 

2  S  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 

0 . 4400000E+08  S  E  -  Young's  modulus  for  skin 

0.2100000  S  u  -  Poisson's  ratio  for  skin 

0.  S  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec**2/in’-4) 

0.  S  ALPHA  -  coefficient  of  thermal  expansion 

0  S  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 

1  S  NSUR  -  control  for  thickness  input  <0  or  1  or  -1) 

N  S  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

N  S  Do  you  want  to  print  out  distributed  loads  along  meridian’ 

H  S 

H  S  SEGMENT  NUMBER  14  14  14  14  14  14  14  14 

H  S  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

5  S  NMESH  -  number  of  node  points  (5  -  min.,  98  -  max.)(14) 

3  S  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 

H  S  REFERENCE  SURFACE  GEOMETRY  FOLLOWS . . . 

1  S  NSHAPE-  indicator  (1,2  or  4)  for  geometry  of  meridian 

6.024000  S  R1  -  radius  at  beginning  of  segment  (see  p.  66) 

38.60300  S  21  -  global  axial  coordinate  at  beginning  of  segment 

5.427000  S  P.2  -  radius  at  end  of  segment 

39.93700  S  Z2  -  global  axial  coordinate  at  end  of  segment 

H  $  IMPERFECTION  SHAPE  FOLLOWS. . . 

0  S  IMP  -  indicator  for  imperfection  (0-none,  1-some) 

H  S  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 

1  S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 

2  S  NZVALU-  number  of  meridional  callouts  for  ref.  surf. 

2  S  NTYPE  -  control  for  meaning  of  callout  (2-z,  3-r) 

38.60300  S  Z( I )  -  axial  coordinate  of  Ith  callout,  z(  1) 

39.93700  S  Z(I)  -  axial  coordinate  of  Ith  callout,  z(  2) 

0.4340000  S  ZVAL  -  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  1) 

0.2750000  S  ZVAL  -  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  2) 

N  S  Do  you  want  to  print  out  r(s),  r'(s),  etc.  for  this  segment? 

B-1 .  Scale  model  housing  data  for  case  4,  Sheet  9. 
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$  DISCRETE  RING  INPUT  FOLLOWS... 

$  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 
$  K=elastic  foundation  modulus  (e.g.  lb/in’-Sjin  this  seg. 

S  LINE  LOAD  INPUT  FOLLOWS... 

S  LINTYP-  indicator  (0,  1,  2  or  3)  for  type  of  line  loads 
$  DISTRIBUTED  LOAD  INPUT  FOLLOWS. . . 

$  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
S  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 
$  NLTYPE-control  (0,1 ,2, 3)  for  type  of  surface  loading 
S  NPSTAT-  number  of  meridional  callouts  for  surface  loading 
$  NLOAD( 1 ) -indicator  for  meridional  traction  (0-none,  1-some) 

S  NLOAD(2)-indicator  for  circumferential  traction 
$  NLOAD ( 3 )  - indicator  for  normal  pressure  (0-none,  1-some) 

S  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

S  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  2) 

$  NTYPE  -  control  for  meaning  of  loading  callout  (2-z,  3-r) 

S  Z ( I )  -  axial  coordinate  of  Ith  loading  callout,  z(  1) 

S  Z(I)  -  axial  coordinate  of  Ith  loading  callout,  z(  2) 

$  SHELL  WALL  CONSTRUCTION  FOLLOWS ... 

$  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 
$  E  -  Young's  modulus  for  sicin 

S  U  -  Poisson's  ratio  for  shin 

$  SM  -mass  density  of  shin  (e.g.  alum. -.00025  lb~sec»*2/in*»4 ) 
S  ALPHA  -  coefficient  of  thermal  expansion 

$  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 

$  NSUR  -  control  for  thichness  input  (0  or  1  or  -1) 

S  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

S  Do  you  want  to  print  out  distributed  loads  along  meridian? 

S 

S  SEGMENT  NUMBER  15  15  15  15  15  15  15  15 

$  NODAL  POINT  DISTRIBUTION  FOLLOWS... 

$  NMESH  =  number  of  node  points  (5  -  min.;  98  -  max.) (15) 

$  NTYPEH-  control  integer  (1  or  3)  for  nodal  point  spacing 
S  REFERENCE  SURFACE  GEOMETRY  FOLLOWS. . . 

S  nshape-  indicator  (1,2  or  4 )  for  geometry  of  meridian 
S  Rl  -  radius  at  beginning  of  segment  (see  p.  66) 

S  Z1  -  global  axial  coordinate  at  beginning  of  segment 

S  R2  -  radius  at  end  of  segment 

S  Z2  -  global  axial  coordinate  at  end  of  segment 

S  RC  -  radius  from  axis  of  rev.  to  center  of  curvature 

S  ZC  -  axial  coordinate  of  center  of  curvature 

S  SROT-indicator  for  direction  of  increasinc  arc  (-1.  or  -1.) 

S  IMPERFECTION  SHAPE  FOLLOWS ... 

S  IMP  =  indicator  for  imperfection  (0-none,  1-some) 

S  REFERENCE  SURFACE  LOCATION  RELATIVE  TO  WALL 
S  NTYPEZ-  control  (1  or  3)  for  reference  surface  location 
$  NZVALU-  number  of  meridional  callouts  for  ref.  surf. 

$  NTYPE  -  control  for  meaning  of  callout  (2-z,  3-r) 

$  Z(I)  -  axial  coordinate  of  Ith  callout,  z(  1) 

$  Z(I)  -  axial  coordinate  of  Ith  callout,  z(  2) 

S  ZVAL  =  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  1) 

$  ZVAL  -  distance  from  leftmost  surf,  to  ref.  surf.,ZVAL(  2) 

S  Do  you  want  to  print  out  r(s),  r'(s),  etc.  for  this  segment? 
$  DISCRETE  RING  INPUT  FOLLOWS ... 

S  NRINGS-  number  (max-20)  of  discrete  rings  in  this  segment 
S  K-elastic  foundation  modulus  (e.g.  lb/in**3)in  this  seg. 

S  LINE  LOAD  INPUT  FOLLOWS. . . 

S  LINTYP-  indicator  (0,  l,  2  or  3)  for  type  of  line  loads 
S  DISTRIBUTED  LOAD  INPUT  FOLLOWS... 

$  IDISAB-  indicator  (0,  1,  2  or  3)  for  load  set  A  and  B 
S  SURFACE  LOAD  INPUT  FOR  LOAD  SET  "A"  FOLLOWS 
S  NLTYPE-control  (0,1,2. 3)  for  type  of  surface  loading 
S  NPSTAT-  number  of  meridional  callouts  for  surface  loading 
$  NL0AD( 1 )-indicator  for  meridional  traction  (0-none,  1-some) 

$  NL0AD( 2) -indicator  for  circumferential  traction 
S  NLOAB( 3 )-indicator  for  normal  pressure  (0-none,  1-some) 


B-1 .  ScaJe  model  housing  data  for  case  4,  Sheet  10. 


-1  $  PN(i)  *  normal  pressure  (p.74)  at  ith  callout,  PN(  1) 

-1  $  PN(i)  -  normal  pressure  (p.74)  at  ith  callout,  PN(  2) 

2  $  NTYPE  «  control  for  meaning  of  loading  callout  (2-2,  3-r) 

39.93700  $  Z(I)  -  axial  coordinate  of  Ith  loading  callout,  2(  1) 

43.29500  $  Z(I)  -  axial  coordinate  of  Ith  loading  callout,  z(  2) 

H  $  SHELL  WALL  CONSTRUCTION  FOLLOWS... 

2  $  NWALL-index  (1,  2,  4,  5,  6,  7,  8)  for  wall  construction 

0. 4400000E+08  S  E  -  Young's  modulus  for  skin 

0.2100000  S  U  =  Poisson's  ratio  for  skin 

0.  $  SM  -mass  density  of  skin  (e.g.  alum. -.00025  lb-sec** 2/in* *4 ) 

0.  S  ALPHA  -  coefficient  of  thermal  expansion 

0  S  NRS  -  control  (0  or  1)  for  addition  of  smeared  stiffeners 

1  $  NSUR  -  control  for  thickness  input  (0  or  1  or  -1) 

N  S  Do  you  want  to  print  out  the  C(i,j)  at  meridional  stations? 

N  $  Do  you  want  to  print  out  distributed  loads  along  meridian? 

H  $ 

H  $ 

H  S  GLOBAL  DATA  BEGINS... 

0  $  NLAST  -  plot  options  (-1-none,  0-geometry,  l-u,v,w) 

N  $  Are  there  any  regions  for  which  you  want  expanded  plots? 

2  S  NOB  »  starting  number  of  circ.  waves  (buckling  analysis) 

2  $  NMINB  -  minimum  number  of  circ.  waves  (buckling  analysis) 

10  $  NMAXB  -  maximum  number  of  circ.  waves  (buckling  analysis) 

1  $  INCRB  -  increment  in  number  of  circ.  waves  (buckling) 

1  $  NVEC  -  number  of  eigenvalues  for  each  wave  number 

0  $  P  -  pressure  or  surface  traction  multiplier 

1  $  DP  -  pressure  or  surface  traction  multiplier  increment 

0  S  TEMP  -  temperature  rise  multiplier 

0  $  DTEMP  -  temperature  rise  multiplier  increment 

0  $  omega  =  angular  vel.  about  axis  of  revolution  (rad/sec) 

0  S  DOMEGA  -  angular  velocity  increment  (rad/ sec) 

H  $  CONSTRAINT  CONDITIONS  FOLLOW.... 

15  $  How  many  segments  in  the  structure? 

H  S 

H  S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  1  1  1  1 

H  S  POLES  INPUT  FOLLOWS... 

1  S  Number  of  poles  (places  where  r=0)  in  SEGMENT 

I  S  I POLE  =  nodal  point  number  of  pole,  IPOLE(  1) 

H  S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS . . . 

C  S  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  S  JUNCTION  CONDITION  INPUT  FOLLOWS... 

N  S  Is  this  segment  joined  to  any  lower-numbered  segments? 

H  S 

H  S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  2  2  2  2 

H  S  POLES  INPUT  FOLLOWS . . . 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS... 

0  S  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  S  JUNCTION  CONDITION  INPUT  FOLLOWS... 

Y  S  Is  this  segment  joined  to  any  lower-numbered  segments? 

1  $  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

1  $  INODE  -  node  in  current  segment  (ISEG)  of  junction,  IN0DE(  1) 

1  S  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

10  S  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  S  IUSTAR-  axial  displacement  (0-not  slaved,  1-slaved) 

1  $  IVSTAR-  circumferential  displacement  (0-not  slaved,  1-slaved) 

1  $  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

1  S  ICHI  -  meridional  rotation  (0-not  slaved,  1-slaved) 

0.  5  D1  -  radial  component  of  juncture  gap 

0.  $  D2  -  axial  component  of  juncture  gap 

Y  S  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

H  S 

H  S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  3  3  3  3 

H  S  POLES  INPUT  FOLLOWS . . . 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

K  5  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS. . . 


B-1 .  Scale  model  housing  data  for  case  4,  Sheet  11 . 
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$  At  how  many  stations  is  this  segment  constrained  to  ground? 

$  JUNCTION  CONDITION  INPUT  FOLLOWS... 

S  Is  this  segment  joined  to  any  lower-numbered  segments? 

$  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

S  INODE  -  node  in  current  segment  (ISEG)  of  junction,  INODE (  1) 

S  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 
S  JNODE  =  node  in  lowest  segmnt  (JSEG)  of  junction 
S  IUSTAR-  axial  displacement  (0-not  slaved,  l-slaved) 

S  IVSTAR-  circumferential  displacement  (0=not  slaved,  l-slaved) 

S  IWSTAR-  radial  displacement  (0-not  slaved,  l-slaved) 

$  I CHI  -  meridional  rotation  (0-not  slaved,  l-slaved) 

$  D1  -  radial  component  of  juncture  gap 
S  D2  =  axial  component  of  juncture  gap 

S  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 
S 

S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  4  4  4  4 

S  POLES  INPUT  FOLLOWS . . . 

S  Number  of  poles  (places  where  r=0)  in  SEGMENT 
$  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS... 

$  At  how  many  stations  is  this  segment  constrained  to  ground? 

$  JUNCTION  CONDITION  INPUT  FOLLOWS... 

$  Is  this  segment  joined  to  any  lower-numbered  segments? 

$  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

$  INODE  -  node  in  current  segment  (ISEG)  of  junction,  IN0D£(  1) 

$  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 
S  JNODE  =  node  in  lowest  segmnt  (JSEG)  of  junction 
C  IUSTAR-  axial  displacement  (0-not  slaved,  l-slaved) 

S  IVSTAR-  circumferential  displacement  (0-not  slaved,  l-slaved) 

S  IWSTAR-  radial  displacement  (0-not  slaved,  l-slaved) 

S  I CHI  =  meridional  rotation  (0-not  slaved,  l-slaved) 

$  D1  -  radial  component  of  juncture  gap 

$  D2  -  axial  component  of  juncture  gap 

$  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 
S 

$  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  5  5  5  5 

5  POLES  INPUT  FOLLOWS. . . 

$  Number  of  poles  (places  where  r-0)  in  SEGMENT 
$  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS . . . 

5  At  how  many  stations  is  this  segment  constrained  to  ground? 

S  JUNCTION  CONDITION  INPUT  FOLLOWS... 

S  Is  this  segment  joined  to  any  lower-numbered  segments? 

$  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

S  INODE  -  node  in  current  segment  (ISEG)  of  junction,  IN0DE(  1) 

$  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 
S  JNODE  =  node  in  lowest  segmnt  (JSEG)  of  junction 
S  IUSTAR-  axial  displacement  (0-not  slaved,  l-slaved) 

$  IVSTAR-  circumferential  displacement  (0-not  slaved,  l-slaved) 

S  IWSTAR-  radial  displacement  (0-not  slaved,  l-slaved) 

$  I CHI  -  meridional  rotation  (0-not  slaved,  l-slaved) 

02  S  D1  =  radial  component  of  juncture  gap 

$  D2  =  axial  component  of  juncture  gap 

S  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  6  6  6  6 

$  POLES  INPUT  FOLLOWS . . . 

S  Number  of  poles  (places  where  r-0)  in  SEGMENT 
S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS . . . 

$  At  how  many  stations  is  this  segment  constrained  to  ground? 

S  JUNCTION  CONDITION  INPUT  FOLLOWS... 

$  Is  this  segment  joined  to  any  lower-numbered  segments? 

$  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

S  INODE  -  node  in  current  segment  (ISEG)  of  junction,  IN0DE(  1) 

$  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

$  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

S  IUSTAR-  axial  displacement  (0-not  slaved,  l-slaved) 

$  IVSTAR-  circumferential  displacement  (C-not  slaved,  l-slaved) 


B-1 .  Scale  model  housing  data  for  case  4,  Sheet  12. 
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1  $  IWSTAR-  radial  displacement  (O-not  slaved,  1-slaved) 

1  $  ICHI  -  meridional  rotation  (O-not  slaved,  1-slaved) 

0.2940000  $  D1  -  radial  component  of  juncture  gap 

0.  $  D2  -  axial  component  of  juncture  gap 

Y  $  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

H  $ 

H  S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  7  7  7  7 

H  $  POLES  INPUT  FOLLOWS . . . 

0  $  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS. . . 

0  S  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  $  JUNCTION  CONDITION  INPUT  FOLLOWS... 

Y  ?  Is  this  segment  joined  to  any  lower-numbered  segments? 

1  S  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

1  $  INODE  «  node  in  current  segment  (ISEG)  of  junction,  INODE(  1) 

6  $  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

7  $  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  $  IUsTAR-  axial  displacement  (0-not  slaved,  1-slaved) 

1  $  IVSTAR—  circumferential  displacement  (O-not  slaved,  1-slaved ) 

1  $  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

1  $  ICHI  -  meridional  rotation  (0-not  slaved,  1-slaved) 

0.6000000E-01  $  D1  *  radial  component  of  juncture  gap 

0.  S  D2  -  axial  component  of  juncture  gap 

Y  $  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

H  $ 

H  $  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  8  8  8  8 

H  $  POLES  INPUT  FOLLOWS . . . 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  $  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS. . . 

0  $  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  $  JUNCTION  CONDITION  INPUT  FOLLOWS... 

Y  $  Is  this  segment  joined  to  any  lower-numbered  segments? 

1  $  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

7  $  INODE  -  node  in  current  segment  (ISEG)  of  junction,  IN0DE(  1) 

7  $  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

9  $  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  $  IUSTAR-  axial  displacement  (0-not  slaved,  1-slaved) 

1  $  IVSTAR-  circumferential  displacement  (0-not  slaved,  1-slaved) 

1  S  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

1  S  ICHI  »  meridional  rotation  (0-not  slaved,  1-slaved) 

0.1950000  $  D1  -  radial  component  of  juncture  gap 

0.  $  D2  -  axial  component  of  juncture  gap 

Y  $  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

H  $ 

H  S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  9  9  9  9 

H  $  POLES  INPUT  FOLLOWS . . . 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS. . . 

0  5  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  $  JUNCTION  CONDITION  INPUT  FOLLOWS... 

N  $  Is  this  segment  joined  to  any  lower-numbered  segments? 

H  $ 

H  $  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  10  10  10  10 

H  $  POLES  INPUT  FOLLOWS . . . 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  $  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS. . . 

0  S  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  $  JUNCTION  CONDITION  INPUT  FOLLOWS. . . 

Y  $  Is  this  segment  joined  to  any  lower-numbered  segments? 

2  S  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

1  $  INODE  -  node  in  current  segment  (ISEG)  of  junction,  IN0DE(  1) 

9  $  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

5  $  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  S  IUSTAR-  axial  displacement  (0-not  slaved,  1-slaved) 

1  S  IVSTAR-  circumferential  displacement  (0-not  slaved,  1-slaved) 

1  S  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

B-1 .  Scale  model  housing  data  tor  case  4,  Sheet  13. 
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1  $  ICHI  -  meridional  rotation  (0-not  slaved,  1-slaved) 

0. 6300000E-01  $  D1  -  radial  component  of  juncture  gap 

0.  S  D2  -  axial  component  of  juncture  gap 

Y  $  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

7  $  INODE  -  node  in  current  segment  (ISEG)  of  junction,  INODE(  2) 

5  5  JSEG  *  segment  no.  of  lowest  segment  involved  in  junction 

7  5  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  $  IUSTAR*  axial  displacement  (0-not  slaved,  1-slaved) 

1  $  IVSTAR-  circumferential  displacement  (0-not  slaved,  1-slaved) 

1  $  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

1  $  ICHI  -  meridional  rotation  (0=not  slaved,  1-slaved) 

0.3470000  S  D1  -  radial  component  of  juncture  gap 

C.  S  D2  -  axial  component  of  juncture  gap 

Y  S  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

H  S 

H  S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  11  11  11  11 

H  S  POLES  INPUT  FOLLOWS . . . 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS. . . 

0  $  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  S  JUNCTION  CONDITION  INPUT  FOLLOWS... 

Y  S  Is  this  segment  joined  to  any  lower-numbered  segments? 

1  S  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

1  $  INODE  =  node  in  current  segment  (ISEG)  of  junction,  INODE (  1) 

10  S  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

13  $  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  S  IUSTAR-  axial  displacement  (0-not  slaved,  1-slaved) 

1  $  IVSTAR-  circumferential  displacement  (0-not  slaved,  1-slaved) 

1  $  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

1  S  ICHI  -  meridional  rotation  (0-not  slaved,  1-slaved) 

-0. 6300000E-01  S  D1  -  radial  component  of  juncture  gap 

0.  S  D2  -  axial  component  of  juncture  gap 

Y  S  Is  this  constraint  the  same  for  both  prebuckling  and  buckling’ 

H  S 

H  S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  12  12  12  12 

H  $  POLES  INPUT  FOLLOWS... 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS. . . 

0  S  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  S  JUNCTION  CONDITION  INPUT  FOLLOWS... 

Y  S  Is  this  segment  joined  to  any  lower-numbered  segments? 

1  S  At  how  may  stations  is  this  segment  joined  to  previous  segs.’ 

1  S  INODE  -  node  in  current  segment  (ISEG)  of  junction,  INODE (  1) 

5  S  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

13  S  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  S  IUSTAR-  axial  displacement  (0-not  slaved,  1-slaved) 

1  S  IVSTAR-  circumferential  displacement  (0-not  slaved,  1-slaved) 

1  $  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

1  S  ICHI  -  meridional  rotation  (0-not  slaved,  1-slaved; 

0. 3000000E-02  S  D1  -  radial  component  of  juncture  gap 

0.  S  D2  -  axial  component  of  juncture  gap 

Y  $  Is  this  constraint  the  same  for  both  prebuckling  and  buckling’ 

H  S 

H  $  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  13  13  13  13 

H  S  POLES  INPUT  FOLLOWS... 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS. . . 

0  S  At  how  many  stations  is  this  segment  constrained  to  ground? 

H  S  JUNCTION  CONDITION  INPUT  FOLLOWS... 

Y  S  Is  this  segment  joined  to  any  lower-numbered  segments? 

1  S  At  how  may  stations  is  this  segment  joined  to  previous  segs.’ 

1  S  INODE  -  node  in  current  segment  (ISEG)  of  junction,  IN0DE(  1) 

12  S  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

10  $  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  $  IUSTAP.-  axial  displacement  (0-not  slaved,  1-slaved) 

1  $  IVSTAR-  circumferential  displacement  (0-not  slaved,  1-slaved) 


B-1 .  Scale  model  housing  data  for  case  4,  Sheet  1 4. 
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1  $  IWSTAR-  radial  displacement  (O-not  slaved,  1-slaved) 

0  $  I CHI  -  meridional  rotation  (O-not  slaved,  1-slaved) 

0.2170000  $  D1  -  radial  component  of  juncture  gap 

0.  $  D?  -  social  component  of  juncture  gap 

Y  $  Is  this  constraint  the  ssune  for  both  prebuckling  and  buckling? 

H  S 

H  $  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  14  14  14  14 

H  S  POLES  INPUT  FOLLOWS . . . 

0  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

H  $  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS... 

0  $  At  how  many  stations  is  this  segment  constraired  to  ground? 

H  $  JUNCTION  CONDITION  INPUT  FOLLOWS... 

Y  $  Is  this  segment  joined  to  any  lower-numbered  segments? 

1  $  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

1  $  INODE  -  node  in  current  segment  (ISEG)  of  junction,  IN0DE(  1) 

13  $  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

5  $  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  $  IUSTAR-  axial  displacement  (0-not  slaved,  1-slaved) 

1  $  IVSTAR—  circumferential  displacement  (0-not  slaved,  1-slaved) 

1  S  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

1  $  ICHI  -  meridional  rotation  (0-not  slaved,  1-slaved) 

0.  S  D1  -  radial  component  of  juncture  gap 

0.  $  D2  -  axial  component  of  juncture  gap 

Y  S  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

H  S 

H  S  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  15  15  15  15 

H  $  POLES  INPUT  FOLLOWS .  .  . 

1  S  Number  of  poles  (places  where  r-0)  in  SEGMENT 

10  $  I POLE  -  nodal  point  number  of  pole,  IPOLE(  1) 

H  S  INPUT  FOR  CONSTRAINTS  TO  GROUND  FOLLOWS . . . 

0  $  At  how  many  stations  is  this  segment  constrained  to  ground7 

H  $  JUNCTION  CONDITION  INPUT  FOLLOWS... 

Y  $  Is  this  segment  joined  to  any  lower -numbered  segments7 

1  $  At  how  may  stations  is  this  segment  joined  to  previous  segs.? 

1  S  INODE  -  node  in  current  segment  (ISEG)  of  junction,  INODE (  1) 

14  $  JSEG  -  segment  no.  of  lowest  segment  involved  in  junction 

5  S  JNODE  -  node  in  lowest  segmnt  (JSEG)  of  junction 

1  $  IUSTAR-  axial  displacement  (0-not  slaved,  1-slaved) 

1  $  IVSTAR-  circumferential  displacement  (0-not  slaved,  1-slaved) 

1  S  IWSTAR-  radial  displacement  (0-not  slaved,  1-slaved) 

1  S  ICHI  -  meridional  rotation  (0-not  slaved,  1-slaved) 

0.  S  D1  -  radial  component  of  juncture  gap 

0.  S  D2  -  axial  component  of  juncture  gap 

Y  S  Is  this  constraint  the  same  for  both  prebuckling  and  buckling? 

H  $  RIGID  BODY  CONSTRAINT  INPUT  FOLLOWS  ... 

N  $  Given  existing  constraints,  are  rigid  body  modes  possible? 

H  $  "GLOBAL!"  QUESTIONS  (AT  END  OF  CASE)... 

Y  S  Do  you  want  to  list  output  for  segment (  1) 

Y  S  Do  you  want  to  list  output  for  segment;  2) 

Y  S  Do  you  want  to  list  output  foi  segment;  3) 

Y  S  Do  you  want  to  list  output  for  segment (  4) 

Y  $  Do  you  want  to  list  output  for  segment (  5) 

Y  $  Do  you  want  to  list  output  for  segment (  6) 

Y  $  Do  you  want  to  list  output  for  segment (  7) 

Y  $  Do  you  want  to  list  output  for  segment;  8) 

Y  $  Do  you  want  to  list  output  for  segment (  9) 

Y  S  Do  you  want  to  list  output  for  segment (10) 

Y  $  Do  you  want  to  list  output  for  segment (11) 

Y  $  Do  you  want  to  list  output  for  segment (12) 

Y  $  Do  you  want  to  list  output  for  segment(13) 

Y  S  Do  you  want  to  list  output  for  segment (14) 

Y  S  Do  you  want  to  list  output  for  segment (15) 

N  S  Do  you  want  to  list  forces  in  the  discrete  rings,  if  any? 

B-1 .  Scale  model  housing  data  for  case  4,  Sheet  15. 
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